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Abstract 

O . 

Gravitational waves from the inspiral of a stellar-size black hole to a supermassive black hole 



C3 



can be accurately approximated by a point particle moving in a Kerr background. This paper 
presents progress on finding the electromagnetic and gravitational field of a point particle in a 



. black-hole spacetime and on computing the self-force in a "radiation gauge." The gauge is chosen 

<n : 

to allow one to compute the perturbed metric from a gauge-invariant component ib$ (or tp^) of 

m 

the Weyl tensor and follows earlier work by Chrzanowski and Cohen and Kegeles (we correct a 
minor, but propagating, error in the Cohen-Kegeles formalism). The electromagnetic field tensor 



and vector potential of a static point charge and the perturbed gravitational field of a static point 
mass in a Schwarzschild geometry are found, surprisingly, to have closed-form expressions. The 



o 
o 

gravitational field of a static point charge in the Schwarzschild background must have a strut, but 

5— i , i/jq and ib 4 are smooth except at the particle, and one can find local radiation gauges for which 

W). 

the corresponding spin ±2 parts of the perturbed metric are smooth. Finally a method for finding 



^ ■ the renormalized self-force from the Teukolsky equation is presented. The method is related to the 

Mino, Sasaki, Tanaka and Quinn and Wald (MiSaTaQuWa) renormalization and to the Detweiler- 
Whiting construction of the singular field. It relies on the fact that the renormalized ip$ (or tb^) is 
a sourcefree solution to the Teukolsky equation; and one can therefore reconstruct a nonsingular 
renormalized metric in a radiation gauge. 
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I. INTRODUCTION 



Among the primary targets of the proposed space-based gravitational wave observatory, 
LISA, are waves from stellar black holes spiraling in to supermassive black holes in galactic 
centers. Because the ratio m/M is small, the orbits and gravitational waves from these 
binary systems can be described to high accuracy by a perturbative expansion with m/M as 
the small parameter. In addition, because tidal forces are small, one can model the system as 
a point mass orbiting a black hole. To zeroth orbit in m/M, the orbit is simply a geodesic of 
the black-hole spacetime. To first order, the particle feels a self-force, whose dissipative part 
is the radiation-reaction force that drives the inspiral. The self-force also has a conservative 



a. 



part that alters the phase of the orbit and of the emitted radiation 

Because the dissipative part of the self-force is antisymmetric under the change from 
ingoing to outgoing radiation - from advanced to retarded fields, it can be computed from the 
half-retarded — half-advanced Green's function. Because this Green's function is sourcefree, 
it is regular at the particle. Approximating the self-force by its dissipative part is an adiabatic 
approximation 3|J4, 5j , and several computations of orbits and waveforms have recently been 
carried out jf], Q, U 9| • 

Including the conservative part of the self force is a more difficult problem, because it 
arises from a field (the half-retarded + half-advanced part of the field) that is singular at 
the particle. One must renormalize the perturbed metric near the particle by subtracting 
off its singular part, a field singular at the position of the particle that does not itself 
contribute to the self force. The MiSaTaQuWa prescription for this subtraction, given by 
Mino, Sasaki and Tanaka[lo| and subsequently in a particularly clear form by Quinn and 



Wald 



111 ], is well understood in a Lorenz (Hilbert, deDonder, Lorentz, harmonic) gauge. A 



Lorenz gauge, however, is not well-adapted to the Kerr geometry: Instead of the decoupled, 

separable Teukolsky equation that simplifies black-hole perturbation theory, one must solve 

a system of ten coupled partial differential equations. 

We report here the beginning of a program to compute the self-force in a gauge appro- 

l I 

priate to the separable wave equation, using a formalism due to Chrzanowski 112] and to 

nn 

Cohen and Kegeles [13, [14| (henceforth CCK) to reconstruct the metric, in what is termed 



a radiation gauge, from either of the gauge-invariant components ipo or ?/> 4 that satisfy the 
Teukolsky equation flol. [l^ |. Cohen and Kegeles clarified and made minor corrections to 
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Chrzanowski's work, 1 and Wald gave a more concise derivation Subsequent work on 
inverting the differential equations that give the Hertz potential is reported in [l8| and 19], 



and a first explicit vacuum reconstruction of the metric in a radiation gauge is given by 



Yunes and Gonzalez 20J. In the present paper, we give a first explicit reconstruction of the 
metric of a point particle in a radiation gauge. And we outline how one can use a version 
of the MiSaTaQuWa renormalization due to Detweiler and Whiting [21] to renormalize the 
gauge invariant component ?/>o (or ip±) of the perturbed Weyl tensor and then to reconstruct 
the renormalized metric in a radiation gauge from the renormalized Weyl tensor. (A differ- 
ent renormalization procedure, also based on a radiation gauge, is given by Barack and Ori 
221.) 



Cohen and Kegeles show that the reconstruction of the vector potential from the spin- 
weight ±1 component 0o or 02 of the electromagnetic field tensor F a/3 is closely analogous 
to the reconstruction of the perturbed metric from ip or and we begin with this simpler 
example. After a mathematical introduction, we obtain in Sect. II I II closed form expressions 
for all components (0o, 4>i and 02) of the electromagnetic field tensor of a static charge in 



a Schwarzschild background, a problem initially solved by Copson|23l|. From 02, we obtain 
the vector potential in a radiation gauge. A radiation gauge exists only where there are no 
sources. We observe that, in a region R, in order to obtain from a Hertz potential a smooth 
vector potential in a radiation gauge, spheres in R enclose no charge; we conjecture that this 
is generally true. One can obtain a smooth vector potential in a radiation gauge by adding 
an I = part of the field in another gauge. (And, at least in a Schwarzschild background, 
one can alternatively add an I = field in a radiation gauge that does not arise from a Hertz 
potential.) 

We next (in Sect. IIVI) obtained a closed-form expression for the components tpo and 
■04 of the perturbed Weyl tensor of a static point mass in a Schwarzschild background. 
There is no consistent solution to the perturbed Einstein equation whose source is a static 
point mass: The mass must be supported. The support, however, can be a strut that does 



1 Eq. (1.7) inRef. appears incorrectly to identify the Hertz potential with a component of the perturbed 
Weyl tensor, but the identification is made correctly later in the paper, apart from a missing complex 
conjugation; Eqs. (1.3)-(1.5) seem incorrectly to imply that radiation gauges can be used when matter is 
present. The equations are correctly given by Cohen and Kegeles, apart from the factor-of-two correction 
we make here. 
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not contribute to the spin-weight ±2 components of the Weyl tensor. In reconstructing 
the metric, we again find that the perturbed metric in a radiation gauge that arises from 
a Hertz potential is smooth only if spheres enclose no perturbed mass. Again one can 
obtain a smooth perturbed metric in a radiation gauge by adding an I = part of the 
field in another gauge. This is part of the freedom one has to add an algebraically special 
perturbation without changing the spin-weight ±2 part of the field on which the radiation 
gauge is based. Regardless of the choice of gauge, one cannot simultaneously make the 
metric smooth everywhere inside and outside the radius of the point mass. To obtain a 
smooth metric in a radiation gauge for the spin ±2 part of the metric, one can add a strut 
by adding an algebraically special perturbation - in this case a perturbed C metric having 
a nonzero deficit angle along an axis through the particle. 

The strut is a feature of our static example that does not appear in an inspiral problem, 
and we conclude our sample reconstructions by considering a point mass in flat space, in 
which we choose a null tetrad that is again not centered at the perturbing mass. The 
formalism is now very close to that of the point charge in a Schwarzschild background, in 
which the only singularity in the perturbed metric outside the particle is pure gauge and 
can be removed by adding an I = perturbation in a different gauge. 

Finally, in Sect. Wl\ we note that the form given by Detweiler and Whiting for the singular 
part of the perturbed metric can be used to find the singular part of tp (or ?/> 4 ). The 
renormalized Weyl tensor component, ip ren = ip^ — ipQ ng , satisfies the sourcefree Teukolsky 
equation. As a result, one can find a nonsingular renormalized metric perturbation in a 
radiation gauge. The motion of a particle is then given to first order in m/M by the 
requirement that it move on a geodesic of the renormalized metric. 



II. MATHEMATICAL PRELIMINARIES 



Greek letters early in the alphabet a, (3, . . . will be abstract spacetime indices; letters 
/i, z/, . . . will be concrete indices, labeling components along the tetrad defined in Eq. (J2]) 

below. We adopt the H signature of Newman and Penrose (NP) 

Schwarzschild metric in the form 

„2 



24 



25], writing the 



ds 2 



4dt 2 -Crfr 2 -rW-r 2 sin 2 



where A = r — 2Mr. 
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We primarily use the Kinnersley tetrad 26[], {e M } (numbered following NP notation) 



r 2 1 1 A 

e? = P = — t a + r Q e? = n a = -t a -r a 

1 A 2 2 2r 2 

e° Em" = -={§" + i<j) a ) e2=m a , (2) 
v2 

where we denote by t a and r a the vectors and d r and by # a and <p a the unit vectors 

— and -cL. The derivative operators associated with l a ,n a , m a and m a are, as usual 

r r sin # 

denoted by D, A, 5 and 5, respectively, but with a boldface A to distinguish this symbol 

from A = r 2 — 2Mr. 

In terms of the nonzero spin coefficients of the Kinnersley tetrad, 

1 cote M 1 A 

P = -r P = - a = ^7Tr 1 = ^ " =_ 2^' (3) 

the corresponding nonzero Christoffel symbols have the form 

1*12 = 27 r 2 22 = -2 7 r 3 33 = 2/? T\, = ^ 
r 3 13 = -p r 3 23 = p ^ = ^1 r 2 43 = - P 

r 4 14 = - P r 4 24 = p r 2 34 = - P r 4 43 = -2/? 



r 3 34 = -2(3 r 4 44 = 2/3. 

Here, for example, V\ 2 e ^ = e 2 V^ei". 



= -F^^V-m^) (5) 



The electromagnetic field F a/3 has independent complex components, 

O = F a(] r m P (4) 

1 

2" 

02 = Faprffn 13 , (6) 

with spin-weights 1, 0, and -1, respectively, where each occurrence of m a (fn a ) contributes 
1 (-1) to the spin- weight. Similarly, the Weyl tensor C a ^s has independent components 



^0 = 


-C^PmPVm 6 


*1 = 




^2 = 




$ 3 = 


—Ca^sPn^rrPn 6 


\J/ 4 = 


—Capys^m^rPm 6 
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(7) 



with spin-weights 2, 1, 0, —1 and —2, respectively. 

In the Schwarzschild geometry, only ^2 is nonzero, and it has value 

M 

^2 = r- 



(8) 



We define the perturbation h a a of a background metric g a a by considering a family of 



metrics g a /3((), and writing 



h 



af3 



A. 



9af3 



(9) 



C=o 



The corresponding components of the perturbed Weyl tensor along the unperturbed tetrad 
will be denoted by the lower-case symbols, i/; , . . . , t])^. For example, 



■03 — — -^Ca/3-yS 



m'n s . 



(10) 



C=o 



Tensor components of the Kinnersley tetrad with spin-weight s have angular behavior 
given by spin-weighted spherical harmonics, s Yi m (6, <j)), when the tensor belongs to an l,m 
representation of the rotation group [27]. To define s Yi m} one first introduces operators 
9 (edth) and 5 that respectively raise and lower by 1 the spin-weight of a quantity r\ of 
spin- weight s: 



drj = — (sin 9) s (dg + i esc 98^) (sin 6) s r/ 
= — (dg + i csc 63^ — s cot 6) rj, 



(11) 
(12) 



drj = — (sin^) s (dg — i csc 98^) (sin 9) s r] 
= — (dg — i csc 9ds + s cot 9) rj. 



(13) 
(14) 



Then, for each value of s, the spin-weighted spherical harmonics are a complete set of 
orthonormal functions on the two- sphere, given by 



s 1 Irn 



[(l-s)\/(l + s)\] 1/2 d s Y lm , 0<s<l, 
(_!)- [(/ + s )\/(i - s )!f/2 Q- Ylm , -l<s<0. 



(15) 
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They satisfy the identities, 



•y* 
s 1 Im 



9 sYlm 



' i \m+s y 

v 1 J —s 1 Imi 



[(l-s)(l + s + l)] 1/2 s+1 Y lm , 
-[(l + s)(l-s + l)] 1/2 s ^Y lm , 
-(I - s)(l + 8 + 1) s Yi mj 

<5(cos0'- cos 0)5(0' -0), 
6u'6 mm >. 



E E * 

Z=|s| m=-l 



/ml 17 ) ' 



S Yl m (9, ' 



(16) 
(17) 
(18) 
(19) 

(20) 
(21) 



The spin ±1 components 0o an d 02 of the electromagnetic field and the spin ±2 com- 
ponents ipo and t/> 4 of the perturbed Weyl tensor satisfy decoupled wave equations, namely 
the Bardeen-Press equation 



and its electromagnetic analog. These are the a=0, spin ±1 



and spin ±2 cases of the Teukolsky equation 



m, 



with 



L 



, M 1 . 
2s ( A "r ]9t + L 



16], 2 and they have the form 
ip = 4ttT, 







dr 



A 



dr 



r 



- 2 m, 



(22) 
(23) 



where ip is any of the functions listed in the first column of Tableland T is the corresponding 
source listed in the third column. The source involves components of the current four-vector 
J a for electromagnetism and of the stress-energy tensor T a/3 for gravity. 

TABLE I: Source function 





S 


T 


00 


1 


<5Ji — (D — 3p) J 3 




-1 


p- 2 {(A + 3p)J±-5J 2 } 




2 


2(6-2/3)[(D-2p)T 13 -ST n ]+ 
2(D - 5p)[(8 - 2/?)T 13 - (D — p)T 33 ] 




-2 


2p" 4 (A + 2 7 + 5p)[(5 - 2p)T 24 - (A + ^)T 44 ] + 
2p~ i (5 - 2/3) [A + 2 7 + 2^)T 24 - 5T 22 ] 



Teukolsky's expressions for the s = ±2 Schwarzschild source functions in Ref. 15j differ from those in 



Ref. 



161 ] by an overall sign. With our conventions, the signs agree with Ref. 16 1 
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Static Green's function 

To compute the fields of static sources, we use the time-independent forms of these 
equations, 

hip = 4ttT. (24) 

Let x and x' denote points of a t =const hypersurface of the Schwarzschild geometry. 
We find below that the source terms for a static particle are related by the operator 9 to a 
point source of the form 5 3 (x, x'), normalized by 



1 = j dV'5 3 (x,x') = j d^'VW) S 3 (x,x'). (25) 
In Schwarzschild coordinates, with yj^g = r 3 A _1//2 sin 9, we have 

5 3 {x, x') = r- 3 A 1/2 5{r - r')5{cos6 - cos 6')5(<P - (j)'). (26) 

Denote by G(x, x') the Green's function satisfying Eq. f[2"4"l) with this source, 

LG(x,x') = -Ati5 3 (x,x'). (27) 

A solution ip(x) to Eq. f l24l then has the form 

${x) = J G(x,x')T(x')dV. (28) 

We construct G(x, x') from solutions to the homogenous equation L-0 = 0. From Eq. ([2"3l . 
these have the form R(r)S(6, 0), with 

MS = \S, ^1(& S+1 ™)+±R=0. (29) 



From Eq. (fl9l) . the solutions to the angular equation are the spin- weighted spherical har- 
monics with eigenvalues A = — (/ — s)(l + s + 1). The radial equation then takes the form 

A " 9 ( A s+1 ^f] -hl-s)(l + s + 1)R = 0. (30) 



Its solutions are given in terms of the associated Legendre polynomials by 

«M = zk<s) + ^«(s)- (31) 

where a and b are arbitrary constants and r = r — M. 
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Using the Wronskian 



W[P?{z),Qi(z)] 



:-iy(i+s) 



(l-z 2 )(l-s)V 

for integral values of I and s, we obtain the Wronskian of the two radial solutions, 



(32) 



W 



( t \ 1 



A 5 / 2 ' \M/' A s / 2 



Qi 



M 



-1 



>s+l 



M 2 (l + s )\ 
A s+1 (l-s)V 



(33) 



Using this relation and the completeness (120]) and orthogonality (l2Tj) relations for the s Y; m 
to compute for each l,m the discontinuity in Eq. ff27j) across r = r', we find for the Green's 
function the form 



(34) 



where r< = min(r, r') — M and r> = max(r, r') — M. 



A. Tetrads smooth at = and 9 = tt and a criterion for smoothness of a tensor 
field. 

The Kinnersley tetrad is singular when 8 = and 9 = tt. In order to disentangle (1) 
the singularity arising from the choice of tetrad from (2) a singularity arising from the use 
of the radiation gauge and (3) a physical singularity associated with a static particle on 
a Schwarzschild background, we will use, in addition to the Kinnersley tetrad, a closely 
related tetrad, {e+} = {P,n Q ,m" = e^m a ,m" = e~^m a }, that is smooth everywhere 
except 6 = 7i. (By smooth we always mean C°°). One can similarly replace m a by e~ l ^m a to 
obtain a null tetrad {e^ }, smooth everywhere except 6 = 0. Because there is no continuous 
vector field nonzero everywhere on a two-sphere, and because m a must be tangent to the 
symmetry spheres if it is orthogonal to the principal null directions, no null tetrad based on 
the principal null directions is smooth everywhere. 

After showing that e ±l ^m a is smooth near 9 = (6 = tt), we will find a simple criterion 
for smoothness of a tensor field near 9 = and 9 = tt that involves only its components in 
the Kinnersley tetrad. 
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To see that e'^m" is smooth near 9 = 0, we introduce a chart {t, X, Y, Z} that is smooth 
on the spacetime exterior to the horizon: With 

X = rsin#cos0, Y = r sin ^ sin 0, Z = rcos9, (35) 

the Schwarzschild metric has the form 

ds 2 = (l - dt 2 - dX 2 - dY 2 - dZ 2 - r2(r 2M 2M) ( XdX + YdY + ZdZ ?i 

smooth for r > 2M. Then e l ^m a is smooth, because its components in this Cartesian chart 
are smooth. We have 

m' z = eV = eVdeZ = -i=^±^, (36a) 

V2 r 



ml* = = e *-L (do + -^A ] X = 4= 

m' Y = e»m x = e»-±r (d e + -^d,) Y = ± 



l-X(X + iY)- r - Z 



r J X 2 + Y 2 

i -i.v • ;y)- ' " z 



r ' X 2 + Y 2 



(36b) 
(36c) 



Finally, ^ 2 _[_ y 2 * S smo °^' 1 ever Y w here except the negative Z-axis, because it is analytic in 
{X, Y, Z}: This is obvious for (X, Y) + 0. To show that it is true for (X, Y) — 0, Z > 0, 
write w = (X/r) 2 + (Y/r) 2 . Then, for Z > 0, 

c,/ n r-Z ^ n i (i _ I) (2 - I) (n-i) 

(37) 

implying F(w) analytic for \w\ < 1. Thus {l a , n a , m® , } is a smooth (and analytic) basis 
for r > 2M, except on the negative Z-axis. Similarly, {l a , n a , m°_ , m"} is a smooth (and 
analytic) basis for r > 2M, except on the positive Z-axis. 

A tensor field is smooth if and only if its components in a smooth basis are smooth. 
Thus a tensor field is smooth on the positive (negative) Z-axis if and only if its components 
along the basis {e^} are smooth. Now a function fe lTn ^ (with / independent of 0) is 
smooth at 6 = (and fe~ irn ^ is smooth at 6 = 7r) if and only if / = g sin" 1 6, with g 
smooth. If / is a spin-weight m component of a tensor in the Kinnersley basis, then je ±tm ^ 
is the corresponding component in the basis {e^}. Replacing sin" 1 9 by 9 m for 9 near and 
(ir — 9) m for 9 near n, we obtain the following criterion for the smoothness of a tensor: 
Proposition. A tensor is smooth near 9 = if and only if its components with spin-weight 
m in a Kinnersley tetrad have the form 9 m g, with g smooth near 9 — 0] and it is smooth 
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near 9 = tt if and only if its components in a Kinnersley tetrad have the form (tt — 0) m g, 
with g a function smooth near 9 = tt. 



In the next three sections, we will use the CCK formalism to write the vector potential 
A a in terms of a scalar $ and to write the perturbed metric h a p in terms of a scalar a 
Hertz potential. We show in Appendix [A] that the Hertz potentials \I/ and $ associated with 
a Kinnersley tetrad are related by phases, e~ 2t ^ and e~ 1 ^ respectively, to the Hertz potential 
associated with a smooth basis. It follows that h a p and A a are smooth if sin 2 9 and 
$/ sin# are smooth. (As noted below, the converse is not true: A singular Hertz potential 
can yield a smooth h or A.) 



III. STATIC CHARGE IN A SCHWARZSCHILD SPACETIME 



A. Outline 



We can now compute the electromagnetic field of a static point charge in a Schwarzschild 
background, finding its complex components (jHEJ) along the Kinnersley tetrad. We first use 
the decoupled equation fl22|) to find fa, obtaining a closed- form expression. Then, following 



Cohen and Kegeles 13, Il4j . we construct a potential in terms of which we compute a 
vector potential A a and the remaining components fa and 0i of the field. 3 Because fa has 
spin- weight -1, it has no / = part, and the monopole part of the field - a change in the 
charge of the black hole - must be determined separately. In this formalism, the potential 
<3? is obtained from fa by solving a second-order differential equation, and freedom to add 
electric and magnetic charge to the black hole can be identified with one of the constants 
of integration. Because the gauges for A a introduced by Cohen and Kegeles exist globablly 
only for a source-free solution, we find local gauges that are singular on different parts of the 
axis of symmetry. Explicit gauge transformations are found relating the vector potential to 
that found by Copson in 1928 [23(] . 



3 Cohen and Kegeles denote by the same symbol ^ Hertz potentials both for the electromagnetic field 
and for the perturbed Weyl tensor. To avoid this ambiguity, we denote by 4> the electromagnetic Hertz 
potential, retaining \P for the gravitational Hertz potential: <£> = ^'(Cohen-Kegeles, electromagnetic), 
^ = ^(Cohen-Kegeles, gravitational). 
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B. Computing 02 from the decoupled wave equation 



Consider a static charge e at a point with radial coordinate 7*0, outside the horizon of a 

Schwarzschild black hole; and choose the 8 = axis to pass through the charge. The charge 

r 

has current 4- vector J a = p e u a , with 4- velocity u a = , t a and with charge density given 

A / 

in the notation of Eq. (I25I) by 

A 1 / 2 

p e (x) = e5 3 (x,x ) = e—^5(r - r )5(cos9 - 1)5(0)- (38) 

The component 02 satisfies the wave equation (1221) with spin s = — 1. Note first that 
in the source term corresponding to s = — 1 in Tabled the angular component J4 vanishes. 
The only contribution to T is then from the component J2 = J a n a , 

T = -\dJ 2 . (39) 
P 

Acting on the spin-weight component J2, the operator 5 is =5, implying 

ry2 

T = -^8j 2 . (40) 



From Eq. (j2J), we have 



r 



1 \ 1AV2 



J2 = ^A^S> 2 r = -2— pe 



±e^5(r-r )5(co S 6- 1)5(0). (41) 



We next decompose the source into spin-weighted spherical harmonics, using the com 
pleteness relation (1201) : 

00 1 

. -5(V — rn) 
2 r 

z=o 

00 i 



1=0 m =-l 



T(x) = --^t-5(r-r )dJ2 £ MM)*£»(0. °) 

Z v Zr 1=0 m=-l 



eA x ' 



2=1 m=—l 



- X Y lm {9, 0)^(0,0). (43) 



It is now straightforward to compute 02 from the Green's function a/) of Eq. 
using orthonormality f|2 lj) of the s 3i m . We have 
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p 2 4>2 



G(x,x')T(x')^jd 3 x' 

A / a \ 1/2 oo I 

4,/M EE 

1=1 m= 

oo n 



X 



M \A 
eA' 



-P 



M 



M 



n 

^^5(r' - r ) £ E + X )] 1/2 ^) W °) 

v n=l p=— n 



i -in 

r' 2 dr'dQ! 



47re(AA 
2 v / 2Mr 



iV2 



oo Z 



S [/(/ + l)]V2 Pi 



lf X < 



M 



Q 



M 



^(9,^(0,0), (44) 



1=1 m=—l 

where r< = min(r, tq) — M, t> = max(r, tq) — M and Ao = r 2 . — 2Mtq. Summing this series 
(details in the Appendix [Bj yields the simple closed-form expression, 

A A sin 9 



2V2r r 2 R 3 



(45) 



where 



R(r, 



= (t 2 + ?l - 2tt cos 9 — M 2 sin 2 9) 1/2 , (46) 
x = r — M and t = r - M. (47) 
An entirely analogous computation of <po, from Eq. (1221) with spin- weight 1, yields the 



expression 



An sin 9 



an sill c , , „. 

""^-ff- (48) 

Our aim, however, is to parallel the Cohen-Kegeles treatment of gravitational perturbations, 
in which the metric and the remaining components of the Weyl tensor are constructed from 
a Hertz potential Here the vector potential A a plays the role of the metric, and A a , 0o 
and 0i are constructed from an analogous Cohen-Kegeles Hertz potential $. 

It is worth pointing out here that, although the electromagnetic field is smooth (C°°) 
outside the particle, its components </>o and 02 are not smooth scalars. 



C. The Cohen-Kegeles formalism for electromagnetism 

We begin with a review of the Cohen-Kegeles Q] formalism for electromagnetic fields 
on type D spacetimes. Underlying the formalism is the following relation between a vector 
potential satisfying the sourcefree Maxwell equation and a scalar $ : 
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Proposition 1. On a vacuum type-D spacetime, let A a be a smooth vector field of the form 

A a = -m a (D + 2t + p)$ + l a {5 + 2(3 + r)$ + c.c, (49) 

where $ satisfies the sourcefree Teukolsky equation for s — — 1 (for 2 ), 

[(A + 7 - 7 + /2)(D + 2e + p) - (5 + a + - f)(8 + 2fi + r)] $ = 0. (50) 

Then A a satisfies the sourcefree Maxwell equations, and the complex components of the 
corresponding field tensor, F a p = VpA a — V a A/3, are given by 

0o = —(D — e + e — p)(D + 2e + (51a) 
20! = -[(D + e + 6 + p-p)(5 + 2p + f) 

+(5-a + p-7r-T)(D + 2e + p)]$, (51b) 

2 = - [(5 + a + P-f)(5 + 2p + f) - \(D + 2e + p)] $. (51c) 

The vector potential given here is in what Chrzanowski and several subsequent authors 
call the ingoing radiation gauge (IRG), in which A a l a = 0. The terminology, however, is 
misleading: An outgoing radiation field in a Lorenz gauge has, near future null infinity, 



wave-vector l a orthogonal to A a (Chrzanowski incorrectly claims in Sect. IV of 12] that 
this is true of an ingoing field near past null infinity). 

An identical proposition holds for the corresponding "ORG" gauge, in which the roles of 
the outgoing and ingoing null vectors, l a and n a , are exchanged. Because of the incorrect 
identification of "outgoing" and "ingoing" in the literature we will simply use the term 
radiation gauges. 

Two comments on what the proposition does not imply: Note first that the component 
02 of the electromagnetic field constructed from A a is not the function $, although both 
$ and 0o satisfy the same sourcefree Teukolsky equation. Second, the fact that 02, say, 
satisfies the sourcefree Teukolsky equation does not guarantee that $ satisfies the sourcefree 
Teukolsky equation. Freedom arising from constants of integration must be used to obtain 
a $ that does so. 

The expressions for A a and for the components of the field tensor in terms of $ are much 
simpler in a Schwarzschild geometry. In this case, all spin coefficients are real and most 
vanish: 

e = r = A = tt = a + (3 = 0. (52) 
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Eqs. ( )49ll51cp then become 

A a = -m a (D + p)<$> + l a (5 + 2(3)<$> + c.c, (53) 
[(A + p)(D + p) - 5(5 + 2(3)] $ = 0, (54) 



-D 2 $, 

-D(6 + 2/3)$, 
-5(5 + 2/3)$. 



(55a) 
(55b) 
(55c) 



For future reference, we note that each operator involving 5 or 5 is equal to 3 or S up to 
factors of —r\[2. In particular, Eq. (I55cl) for 2 can be written in the form, 



2r 2 



(56) 



D. The scalar potential 

We now invert the CK relation (1561) to find $ in terms of the field component <p2 of a static 
charge, using the closed-form expression for <p2- From Eq. fflBl and the axisymmetry of 
2 , we see that the leftmost 8 is just — dg. Integrating with respect to 9, we have 

e A A 



5$ 



V2r 



V2 



r L 



sin 9 
R 3 (r,9) 
M 2 cos 9 - tt 
R 



d6 



+ a(r) 



(57) 



with a(r) an arbitrary function of r. Similarly, after writing 5$ = — (sin#) 1 dg(sin9$), a 
second integration with respect to 9 yields 



y/2 r sin 9 



M 2 cos 9 - tt 
R 



+ a(r) 



dcos9 



[R + a(r) cos 9 + b(r) 



\[2 r sin 9 
with b(r) again arbitrary. 

Eq. fl58|) has the form of a particular + a homogeneous solution to Eq. fl56l) . 



(58) 



e 



\/2 r sin 9 ' 

e a(r) cos 9 + b(r) 
\[2 r sin 9 



(59) 
(60) 

(61) 
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The particular solution $p already satisfies the sourcefree Teukolsky equation. The homoge- 
neous solution, however, does so if and only if a"(r) = b"(r) = 0. With this restriction, 
the solutions constitute a 4-parameter set, specified by 

a(r) = ao + air, b(r) = &o + (62) 

By construction, a and b encode no information about 02. Because we expect 02 to carry all 
information about the I > 1 part of the field, the value of 0o should similarly be independent 
of a and b. This follows directly from Eqs. (I55ap and (!62l . which together imply that 
0o [&h] = 0. Then a and b can carry information only about the monopole field. As we show 
in the next section, the parameter ao corresponds to adding charge to the black hole. The 
remaining parameters ai, bo and b\ correspond simply to gauge transformations (with b\ the 
trivial gauge transformation, altering neither the field nor the vector potential). 

The potential $ is in general singular on the z axis because of an overall factor of 1/ sin 9. 
This singularity is carried to the vector potential as well, because the angular pieces of 
equation (1531) contain only radial derivatives. We will see that one can choose a\ and bo to 
make A a smooth on the axis for r > ro or for r < ro, but not both. 

E. Completion of the solution 

Fields corresponding to &h 

The vector potential Aq,[$#], associated with the homogeneous part of the potential 

$ H = ^ e ■ n t( fl o + fl i r ) cos 6 + ( 6 o + hr)] , (63) 
V 2 r sm 



has in the radiation gauge the value 



i ., ,„cote \ V2 

-ao —l a — v 2 m a — a\l a + o — ; — -p.fn a 



(-L-V2' 

\r 



r I r sm U 



+ c.c. (64) 



Using the relations 



1 1 / 6 
l a = Vq,m, — ; — -m a = pV a I log tan - + icj) ) , (65) 



rsinfl a ^2 V b 2 
with u the null coordinate u = t — r — 2Mln(r/2M — 1), we have 



= - a ow— (fa - v 7 ^ cot#m a ) +c.c. + V Q A, (66) 



2r r 
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where the gauge scalar A is given by 

' " X fee. (67) 



2r 



Q 

a\U — b ( In tan — + i<p 



Thus, as claimed, the vector potentials associated with ai and b are pure gauge, and that 
associated with b\ vanishes. 

The identification of with a change in the black hole's charge (electric and magnetic) 
follows from the form of the corresponding field tensor: By a test charge q on a Schwarzschild 

black hole is meant a radial electric field, with F a pi OL r (i = — . Equivalently, only the spin-0 

r q 

part of the field is nonzero, and it is real, with value <f>\ = -. From Eqs. ()55aH55cl) , 

the field associated with Q H for a real is the field of a black-hole test charge of magnitude 
Q = ^: 

»"0 

M®h] = U^h] = 0, (68a) 

M*h] = -^r-r (68b) 
Ztq r z 

Similarly, a test magnetic charge g on a, Schwarzschild black hole has the meaning of a radial 
magnetic field with F a ^Q a <p^ = Equivalently, only <pi is nonzero, and it is imaginary, 

with value X = — -. And the field associated with for a® imaginary is the field of a test 
2r 2 

magnetic charge of magnitude g = " ' 



If we replace a in by Q = r a /e, we can summarize the last paragraph as follows: 

Q 

The electromagnetic field associated with $ = — j= cot 9 is a monopole field with electric 
charge Re Q and magnetic charge Im (-Q). 



Fields corresponding to $p 

The vector potential corresponding to $p has the value 



, r = n e fx— rncos^ R\ . , e , , , 

= (— 7) {m « +m °> + t^r <"» " cos ) «. 

The corresponding components of the field, given by Eqs. (l55aH55cl) . are 

A sin 6* 

e M 2 cos 6 — rr 



(69) 



(70a) 



*i = H^ — s-^;- (70b) 

A ° A8i "" (70c) 



2V2r„ 
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As we now show, the black-hole charge for this field is nonzero. To find the field of a point 
charge for which the black hole has zero charge, one must add to </>i[$p] a multiple of 
that has black-hole charge opposite to that of our particular solution. The total charge 
of the particular solution (black-hole charge + point charge) can be found by evaluating 
AtiQ = j s E ■ dS = j s 2<p 1 dS over any two-sphere enclosing the black hole and particle. We 
can, for simplicity, take the sphere to be at spatial infinity, writing 

Q = _L to / a,. f^ 2 ^-M r , rfSi 

4ir^«J Sr \2r r R J 

^0 / 

Thus the particular linearized solution associated with $p describes a particle of charge e 
outside a black hole of charge 

M 

Qbh = -e— . (72) 
r 

We will see in the next section that any homogeneous solution that makes the vector 
potential smooth for r < r has a black hole charge that exactly cancels the black hole 
charge of the particular solution: A radiation gauge is smooth for 2M < r < only if 
spheres with 2M < r < r$ have no enclosed charge. Similarly, we will see that a radiation 
gauge is smooth for r > r only if the total charge (black hole + particle) vanishes: Spheres 
with r > ro enclose no charge. 

Parameter choices for solutions smooth on the axis of symmetry 

As we showed in Appendix |AJ a sufficient condition for A a to be smooth is that the 
Hertz potential $ is smooth; and $ of the form given by Eqs. (1591) - (|62|) is smooth if 
and only if $ = 0(6) near 6 = and $ = 0(n — 6) near 6 = 7T. Smoothness of $, 
however, is not a necessary condition for smoothness of A a , because $ involves a parameter 
b± that corresponds to a vanishing vector potential. By omitting conditions on bi, we obtain 
necessary and sufficient conditions for smoothness of A a . (It is slightly more efficient to find 
the conditions from $, but we will also obtain them directly from the components of A a as 
a check.) 

From Eqs. (!60l) and (|63|) . we have 

$ = $ p + $ H = —1 -[R + (a + a x r) cos# + b + b x r\. (73) 

V2 r sind 

Smoothness on the 6 = part of the axis (the part of the axis that contains the particle) is 
enforced by different conditions below and above the particle's position at r = r . Formally, 
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the difference arises from the fact that the quantity R = (r 2 + tg — 2tto cos 6 — M 2 sin 2 #) 1//2 
in the numerator of the above expression for 5> has the value R{9 = 0) = \r — Tq\. Then $ 
is smooth at 6 = if and only if cto + a\r + 6 + b±r = — \r — r Q \, or 

ao + b = -t , Oi + 6i = 1, r < r (74a) 
«o + &o = r o, Oi + fei = -1, r > r . (74b) 

Similarly, from the value R(8 = tt) = r + ro — 2M, it follows that $ is smooth at 6 1 = it if 
and only if 

a -b = r - 2M. (75) 

The conditions on and 6o i n EcjS- (1741) are inconsistent, implying that no choice of param- 
eters yields a radiation gauge smooth both below and above the particle. One can at most 
choose parameters that satisfy condition (1751) and one of the two conditions (174|) . These 
choices yield a vector potential A a 

(I) smooth everywhere except 6 = 0, r > r , or 

(II) smooth everywhere except 8 = 0, r < r . 

The implied choice of parameters for each case is given in Table [III Note that any choice 
TABLE II: Singularity location for two choices of parameters 



Parameter 


Singularity at 


Singularity at 




6 = 0; r > ro 


= 0; r < r 


a 


—M 


r - M 


bo 


-(ro-M) 


M 



of the parameter a\ is permitted, because one can choose the trivial parameter b\ to satisfy 
the second condition without altering the value A a . 

In case (I) the vector potential is, in particular, smooth in the exterior Schwarzschild 
geometry for all r < r . The parameter value = M corresponds to a black hole charge 
eMj r that exactly cancels the black hole charge (172]) of the particular solution. In case (II) 
the vector potential is smooth for all r > tq, and the parameter value ao = tq—M corresponds 
to a black hole charge — e(l — M/r ) that cancels the total charge of the particular solution: 
The full solution has a black hole with charge — e and total charge 0. Thus as noted above, 
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a Hertz potential $ yields a smooth radiation gauge on a sphere only if the total electric 
and magnetic flux through the sphere vanish. 

This is not surprising. The radiation gauges are designed for sourcefree solutions and are 
constructed from the spin ±1 parts of the field - in our case from parts of the field with 
no I = part. For a point-charge with no charge on the black hole, we can find regular 
radiation gauges outside and inside the charge for the I > 1 part of the field and can add to 
them the vector potential in another gauge for the Coulomb I = field. 

It makes sense to use a radiation gauge to describe the part of the field with spin greater 
than and to add an / = part in another gauge. In particular, by adding the solution 
describing a perturbation of the black hole's charge in a Coulomb gauge, namely 

A - 9. 

r 

0/1 r 2 \ 
A a = A t V a t = j ( - l a + —n a J , (76) 

we obtain a solution with arbitrary black hole charge, whose gauge singularity is restricted 
to lie on the 9 = ray, either above or below the particle. 4 

The vector potentials corresponding to cases (I) (singular only for 9 = 0, r > r ) and 
(II) (singular only for 9 = 0,r < r ) are, respectively, 

„ m e / rt - M 2 cos 9 \, 
A { ' ] = — — + M )l 



r r \ R 

I e ( R-^-Mcos9 r-tpcosjA , (77) 

1 7^ 7 I 5 ] \jn a -+- m a ), {<<& 

y/2r sm0 \ r R J 

^ = \ 

r r V R J 

V2r sm9 \ r R J 

Here we have set the parameter a\ to zero. 

As we have noted, Eqs. (I77ap and (177bl) are physically different solutions, with the former 
having an uncharged black hole, the latter a black hole with charge — e. Adding to the vector 
potential f!77bj) the I = (Coulomb) perturbation (176]) that cancels the charge — e on the 



4 Although one cannot obtain the I = vector potential in a radiation gauge by using a Hertz potential of 
the form (|73|) . the vector potential can be expressed in a radiation gauge: Writing / = Qln {^J^J ~ l)> 



Q 

we have A a = A a — V a / = — l a , implying A a l a = 0. 

r 
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black hole, we obtain a vector potential that is gauge related to that of Eq. (I77aj) and is 
smooth except on the axis below the particle: 



= — 

r r 

+ 



e / tt - M 2 cos 9 1 A c 



R 



2 r 



er 



A 



V2 r sin 



i? + M + t cos 6> t — t cos 9 
r R 



(m a + m a ). 



(78) 



Finally, as a check, we directly verify the smoothness conditions governing ao and bo in 
flU]) and (175]) . using the components of A a [$] = A Q [$ P ] + A Q [$ H ] given by Eqs. flMD and 
([6T?]) . Because v4 2 is already smooth for arbitrary values of the parameters, A a will be smooth 
if and only if A4 is 0(9) near = and 0(n — 9) near 9 = tt. Near = 0, we have 

e a + 6 + ^0 



A 4 [$] 



\/2r r 9 
e a + bo-r 

\^2r r 9 



+ 0(9), r<r 
+ 0(9), r>r . 



(79a) 
(79b) 



Near the axis of symmetry on the side of the black hole opposite to the particle (9 = 71) 
the corresponding expansion of A a has the form 

e -Oq + b + r -2M 



A> 



+ 0[( 



71 



(80) 



a/2 r r 7T-9 

The conditions on ao and 60 in Eqs. ( 1741) and (175]) are exactly conditions that the coefficients 
of 9~ x and (tt — 9)~ l vanish. 

We can explicitly verify that the vector potentials A^P and A^ we have obtained for a 
point charge outside an uncharged black hole are related by a gauge transformation to the 
solution obtained by Copson. In particular, using the gauge function 



K(t,r,9, 



— eln 

eM 



t cos 9 — t + R 



+ — In [2 (t cos 9 - t + R)] 

r 



+ 



In 



r 



cos 9 + 1 

A - M 2 sin 2 9 + Mr - t r cos 6» + (M(l + cos 1 



r )R 



M(l + cos9) sin 2 9 



,(81) 



we find A { a ] + V a A = A t V a t, with 



A = — 

rr 



e / tt — M 2 cos 9 



R 



+ M 



(82) 
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IV. STATIC GRAVITATIONAL PERTURBATIONS 



A. Outline 

Notation. Because the components V'o>'01j'03 an d ^4 °f t ne Weyl tensor vanish in the 
background Schwarzschild spacetime, we use the same symbols to denote the components 
of the perturbed Weyl tensor. Because ^2 is nonzero in the background spacetime, we 
denote the perturbed component by ip2- 

Outline of Computation. To compute the linearized gravitational field of a point mass in the 
Cohen-Kegeles formalism, we roughly parallel the calculation of the electromagnetic field of 
a point charge, but there are two primary differences. First, where the electromagnetic field 
tensor is gauge invariant, only the extreme spin components, tpo and ^4, of the perturbed 
Weyl tensor are gauge invariant. Second, the electromagnetic field of a static test charge 
makes sense, but a static test mass must be supported by something that itself contributes 
to the linearized gravitational field, and in our solution this is seen as a conical singularity 
(a deficit angle) along the axis of symmetry inside or outside r = r . Nevertheless, the 
extreme spin components of the perturbed Weyl tensor are smooth everywhere except at 
the position of the particle: The part of the perturbation with \s\ = 2 does not know about 
the strut. As in the electromagnetic case, we obtain simple closed-form expressions for these 
gauge-invariant parts of the static perturbation. 

We first use the decoupled equation (j2"2l to find ip4, again by taking angular derivatives 
of the static Green's function. Using a radiation gauge, we construct a potential in 
terms of which we compute the components h a p of the perturbed metric and the remaining 
components ipi,ip2 and 1P3 of the field. In so doing, we correct of factor-of-two error in the 
Cohen-Kegeles formalism. 

Because ip4 has spin- weight —2, it has no / = or / = 1 parts. The monopole part of 
the field corresponds to a change in the mass or area of the black hole, the dipole part to a 
change in its angular momentum and center of mass. As in the case of the electromagnetic 
charge, the change in the black hole's mass must be specified separately. The change in its 
center of mass appears as a gauge transformation. The potential \1/ is obtained from ^4 by 
solving a fourth-order differential equation, and the freedom to add mass and spin to the 
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black hole can be identified with two of the the constants of integration. We again find local 
gauges that are singular on different parts of the axis of symmetry (the singularities in these 
radiation gauges were identified by Barack and Orij^), with a residual conical singularity 
(the strut) that cannot be removed by a choice of gauge. 



B. Computing ip4 from the decoupled wave equation 

Consider a static point particle of mass m at a point xq with coordinates r = ro, = 0. 
The particle has density 

p m (x) = m5 3 (x,x ) (83) 

and stress-energy tensor 

T aP = p m u a vP, with u a = -^t a . (84) 

In the source term of the Bardeen-Press equation for p~ A ip^ given in Table [H only the 
component T 22 is nonzero, and the source is simply 

T = -2{5- 2(3)5 p~ A T 22 , (85) 

with 

A3/2 

T 22 = T aP n a n? = m— 5(r - r )6(cos6 - 1)5(0). (86) 

By writing the angular derivatives in the source in terms of <3 and using equation (1201) . 
we can expand the source in spin-weighted spherical harmonics: 

T(x) = -2(5- 2(3)5 p~ A T 22 
= -r 2 8 2 T 22 

^3/2 00 1 



— m 



1=0 m=-l 
l= 2 m=-l 



„ 2 Y lm (6, 0)^(0,0). (87) 



.(* -2)!. 

We can now obtain ip^ from the Green's function of Eq. (1341 . Eq. fl28l) for p~ A ip4 has the 
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form, 



P Hi 



G(x,x')T(x')y/^d 3 x' 

oo I 



(88) 



4ttA 



MA 
mA' 



(j ~ 2 ) ! p2 \ n 2 f 1 

(Z + 2)! 1 \m)^ 1 



4r 



1=2 m=-l 

3/2 oo n 

5{r'-r )Y^ Yl 

n=2 p=—n 



/3 



(w + 2)! 
_(n-2)!_ 



1/2 



_ 2 y np (^0')i;;(o,o) 



and, using orthogonality of the spin-weighted harmonics, we have 



1p4 



7rm AA^ 
'~M r 4 r 



1/2 oo I 



E E 

Z=2 m=-; 



(*-2)! 
(Z + 2)! 



1/2 



P 



f (SfW(^) -2^0)^(0,0). (90) 



M 



It is again possible to sum this series to obtain the closed-form expression (details in Ap- 
pendix [B]), 



3Aj /2 A 2 sin 2 6 



(91) 



An entirely analogous computation of ip , from Eq. (122]) with spin- weight 2, yields the 
expression 

3A 3 / 2 ^2tf 

(92) 



3A^ /2 sin 2 9 
ipo = -m — - 



r P 5 



C. The Cohen-Kegeles formalism for gravity 

The formalism of Cohen and Kegeles [14J] again relates ip4 to a scalar potential \I/ which 
we can then use to reconstruct the perturbed metric h a p. In each of the key relations, 
Eqs. dHMHHD below, we correct an error that appears in the Cohen-Kegeles papers and has 
been repeated in the subsequent literature: The right hand side of each equation is reduced 
by the factor 1/2 from the corresponding Cohen-Kegeles expressions. 

Proposition 2. On a vacuum type-D spacetime, let h a p be a smooth tensor field of the form 

Kp = -{Up [(8 + a + 3fi - f)(5 + 40 + 3f) - X(D + 4e + 3p)] 
+m a m p (D - e + 3e - p) (D + 4e + 3p) 
-/ {Q m /3) [(D + e + 3e + p-p)(5 + 4/5 + 3f) + 

(5 - a + 3(3 - 7T - f )(D + 4e + 3p)] + c.c, (93) 
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where \I/ satisfies the sourcefree Teukolsky equation for s = — 2 (for ^4), 

[(5 + 3a + (3 - t)(6 + 4£ + 3f) - (A - 7 + 37 + y)(D + 4e + 3p) + 3^ 2 ] * = 0. (94) 

Then /i a( g satisfies the perturbed vacuum Einstein equation, and the complex components 
of the corresponding perturbed Weyl tensor are given by 



V>o = -{D-3e + e-p){D-2e + 2e-p){D-e + 3e-p){D + 4e + 3p)^, (95) 
V>i = -[(D - e + e + p - p)(D + 2e + p - p)(D + e + 3e + p - p)(5 + + 3f) 

o 

+{D - e + e + p - p) (D + 21 + p - p) {5 - a + 3(3 - n - f ) (D + 4e + 3p) 

+ (D - e + e + p - p){8 - 2a + 2(3 - 2n - f )(D - e + 3e - p)(D + 4e + 3p) 

+ (5 - 3a + (3 - 3tt - f )(D - 2e + 2e - p)(D + e + 3e - p)(£> + 4e + 3p)]f , (96) 

tp 2 = —[{D + e + e + 2p-p){D + 2e + 2e + 2p-p)(5 + a + 3p-f){6 + + 3f) 
+(D + e + e + 2p-p)(5 + 2p-7r-f)(D + e + 3e + p-p)(5 + 4P + 3f) 
+{D + e + e + 2p - p)(5 + 2(3 - n - f )(5 - a + 3(3 - n - f )(D + 4e + 3p) 
+ (S - a + (3 - 2n - f)(5 - 2a + 2(3 - 2n - f)(D - e + 3e - p)(D + 4e + 3p) 
+ (5 - a + (3 - 2tt - f)(D + 2e + p - p)(5 - a + 3(3 - ix - f )(D + 4e + 3p) 
+ (6 - a + - 2tt - f)(D + 2e + p - p)(D + e + 3e + p - p){5 + 4/3 + 3f)]*, (97) 

^3 = -[( J D + 3e + e + 3p-p)(5 + 2a + 2 ( 5-f)(5 + a + 3^-f)(5 + 4^ + 3f) 

o 

+ (5 + a + /3 - 7T - r )(£> + 2e + 2e + 2p - p)(5 + a + 3(3 - f )(<5 + 4(3 + 3f ) 
+ (5 + a + /? - 7T - f )(5 + 2/3 - 7T - f)(D + e + 3e + p - + 4^ + 3f ) 
+ (5 + a + (3 - 7T - f){5 + 2(3-n- f)(S - a + 3(3 - n - f )(D + 4e + 3p)]^, (98) 
i) 4 = -{(5 + 3a + p-f)(5 + 2a + 2p-f){5 + a + 3P-f)(5 + 4p + 3f)f 

+3^ 2 [r(5 + 4a) - p(A + 4 7 ) - p(D + 4e) + tt(5 + 4(3) + 2^ 2 ]^}. (99) 

Note that the metric h a p is tracefree and satisfies h^l 13 = 0. Because these are five real 
conditions on the components of h a p, and one has only four gauge degrees of freedom, one 
cannot fi nd a radiation gauge for a generic metric. As Price, Shankar and Whiting show, 
however 



30, 



311 ] . such a gauge exists locally for vacuum perturbations of any type D vacuum 



spacetime. In particular a perturbation that changes the mass of a Kerr spacetime (with J 
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and M nonzero for the background) cannot be written in a radiation gauge, but radiative 
perturbations and a perturbation that changes the angular momentum of Kerr can be. 

As we show below (and as Price et al. found independently), however, one can express a 
perturbation that changes the mass of a Schwarzschild black hole in a radiation gauge. We 
exhibit two examples acquired from a Hertz potential that are each singular on a ray. Al- 
though the Hertz potential does not yield a form of the mass perturbation that is nonsingular 
in the entire exterior Schwarzschild geometry, the family of radiation gauges includes such 
a perturbation, and we present it below. The Hertz potential formalism, appropriate for 
spin-weight greater than two, yields at least for s = perturbations (for / = perturbations 
of Schwarzschild), a restricted family of radiation gauges. 

By exchanging the roles of the ingoing and outgoing null vectors one obtains analogous 
propositions for a radiation gauge in which h a prfi = 0, h = 0. 

These equations are much simpler on a Schwarzschild background. The form of h a p is 

h a/3 = -[lJ f3 (6 + 2(3)(6 + 4/3) + m a mp (D-p)(D + 3p) 

-2l {a m p) {D + p)(6 + 4/3p + c.c, (100) 

where 

[(6 - 2(3) (5 + 4/3) - ( A + 2 7 + p) (D + 3p) + 3^ 2 ] f = 0. (101) 
The expressions for the complex components of the Weyl tensor become 



V'o = 


2 


(102) 


Vl = 




(103) 


1p2 = 


l -D\5 + 2P){5 + A^ 1 


(104) 


^3 = 


irw(tf + 2/?)(S + 4/?)f, 


(105) 


^4 = 


l -{5- 20)6(6 + 2(3)(5 + 40)# 





--q 2 {pD + p(A + 4-f) -2tf 2 ]tf. (106) 



In the next section, we solve Eq. (I106j) to find the scalar potential ^ from the closed-form 

expression (J9T1) for ip^. In the second line of Eq. ( 11061) . the expression in brackets becomes 

simply ——d t , when one substitutes the Schwarzschild values of the operators and spin 
2r 
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coefficients (J3J), and the value — M/r 3 of the unperturbed ^2- In our static case, the second 
line then vanishes, leaving a set of angular derivatives of Replacing ^4 on the left side 
by its value (l§Tj) . we have 



itaAjpAW* i ( j_ 2/W + 2/ , )( j + 4/3)f (107) 



4r r 4 i? 5 2 



g^*. (108) 



or 



§4 - = _6mAT AW. (109) 
r ix 5 

D. The scalar potential ^ 

We can now begin the integration of Eq. f ll09j) for ty. Because ^4 has spin-weight —2 
and B lowers the spin-weight by 1, ^ has spin-weight +2. From the definition (TTH]) of 5, its 
action on the spin-weight —1 quantity is given by 

B 4 ^ = -sin0d e ( ^-rB 3 ^ ] . (110) 
\sm# / 

Integrating this equation with respect to 8, we have 

6mA 2 A„ /2 „ [ sin 6d6 . . 

= ^ s in# / — 111 

r J R 5 

6mA 2 Af . /■ dcostf 

= — sin0 / rrj (112) 

r o 7 (t 2 + t 2 -2tt cos#-M 2 sin 2 #) 1 

2m . r M 2 cos — rr 



sin g{ C ° S 3 rt ° [r 2 r 2 - 3M 2 (t 2 + t 2 ) + 3M 4 + 4rr M 2 cos - 2M 4 cos 2 9] 
+ a(r)|, (113) 



with a(r) arbitrary. 

Next, acting on the spin- weight quantity 9 2 \1/, 5 has the form 



9 3 # = -d e 3 2 y, (114) 
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and we have 

8 2 * = 2m 



r oA 



/{ 



M 2 cos 9 - tt 



Ay 2 j i & 

+a(r)j<i cos 9 



[th 2 - 3M 2 (r 2 + tg) + 3M 4 + 4rr M 2 cos - 2M 4 cos 2 0] 

(115) 



2m 



r oA 



1/2 



t't 2 , + M 2 (t 2 + tg) — M 4 — 4tt M 2 cos g + 2M 4 cos 2 g 
(t 2 + r 2 - 2tt cos 9 — M 2 sin 2 0) 1/2 

+ a(r) cos# + b(r) 



(116) 



with b(r) arbitrary. 

Continuing in this way, we have 

d 2 ^ = l —d e (sin 9m) , 

sin 9 v 7 

whose integration entails a third free function c(r): 



(117) 



sin» = 



2m 



1/2 



r oA 



2m 



r 2 r 2 + M 2 (r 2 + r 2 ) — M 4 — 4tt M 2 cos + 2M 4 cos 2 9 

R 



+ a(r) cos 9 + b(r) 



dcos9 



The final integration of 



(tt - M 2 cos 9)R + -a(r) cos 2 9 + 6(r) cos 6> + c(r) 



(118) 



sin» = -de (sin 2 ^) , 

sin 9 



yields 



sin 2 9^ 



2m 



2m 



(tt - M 2 cos 6»)i? + -air) cos 2 + 6(r) cos 9 + c(r) 
111 

■-R 3 + -a(r) cos 3 9 + -6(r) cos 2 6> + c(r) cos 6> + d(r) 
3 6 2 



(119) 



dcos9 



(120) 



We thus find a remarkably simple particular solution ^/p, together with a homogeneous 
solution ipn involving four free functions, a(r),b(r),c(r) and d(r): 



\p = vj/p + \S? H 

= 



2 „ R 3 
-A- 



3 sin 2 fl 



1 



sin 2 9 



1 i 

-a(r) cos 3 9 + -6(r) cos 2 + c(r) cos# + d(r) 



(121) 
(122) 

(123) 
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TABLE III: Free functions of the 


gravitational scalar potential 


Function 


Value 


a{r) 


air 2 {r - 3M) + a 2 




b(r) 


bir 2 + (r - M)b 2 


- ^a 2 + c 2 (r - M) + cir 2 


c(r) 


— d(r - M)(r 2 + 4M 2 ) 


d(r) 


1 1 

-&ir 2 + -6 2 ^ + dir 2 (r - 


3M) + d 2 



where we denote by A a constant that appears repeatedly throughout this section, 

Again the particular solution ^ p already satisfies the homogeneous Bardeen-Press equa- 
tion for s — 2. Requiring that also be a solution restricts each free function to a 
polynomial whose form is given in Table IIHI Then the scalars ^ of Eq. (I12ip form an eight- 
complex-parameter set of potentials that generate vacuum perturbations for which ^4 has 
the value (I9T1) . with a%, a^, 61, 62, Ci, C2, di, d 2 the eight complex parameters. 

Now each \l/p generates a vacuum perturbation for which the gauge invariant component 



^4 vanishes. It follows that ipo vanishes as well, and a theorem of Wald [32| restricts the 
perturbations associated with \1/ # to combinations of perturbations that change the mass 
and spin of the black hole, perturb Schwarzschild to a C metric or to Kerr-NUT, or are pure 
gauge. We begin the next section by characterizing in this way the perturbation associated 
with each of the eight parameters of \l/p. 



E. Completion of the solution 

Fields corresponding to 

We will need only real values of the parameters a±, . . . , d 2 , but will also include complex 
values of a 2 , because Im(a2) corresponds to a perturbation taking Schwarzschild to Kerr, 
adding angular momentum to the Schwarzschild geometry. 

We note first that the perturbed metric associated with a± is a C-metric perturbation. 
The components of the perturbed metric and Weyl tensor in the radiation gauge have the 
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following forms: 



h 22 = -aiA2(r - 3M) cos 9, 

AM 2 (r -2M) -r 3 sin 2 9 

h 23 = a x A : : 

y 2 r z sin 

2M (r sin 2 - AM) cos 6 1 

^33 = UlA r-TT, ; 

r sm 9 



(125a) 
(125b) 

(125c) 



^ = 1P2 = ^4 = 0, (126a) 

* - M^, (126b) 

r 4 sm 
, 3Msin# 

^ 3 = -a x A = . (126c) 

2V2 r 2 



To verify that this is a perturbed C metric, we refer in Wald [32] to the perturbed tetrad 
produced by what he calls p. Using this perturbed tetrad to compute the perturbation in the 
metric g a p = 2n^lp) — 2m^fn^, we calculate a perturbed metric. It is then straightforward 
to find a gauge transformation from Wald's gauge to our radiation gauge. 

When a 2 is real, the perturbation associated with it is pure gauge: The nonzero compo- 
nents of the perturbed metric, 

, 2 cos 9 , 
h 22 = -a 2 A—^—, (127a) 

, , \f2 sin6> ,., . 

h 2 3 = a 2 A — , (127b) 

h 33 = 0, (127c) 



and Weyl tensor, 



, 3 sin 9 

V>i = -a 2 A— — , (128a) 



2^2 

are the components of £^g a ^ and L^C^s for the gauge vector 



. 3 cos 9 , 
^2 = a 2 A—^, (128b) 

^3 = a 2 Ap^-, (128c) 



, cos 9 

6 = , (129a) 

, (r + 2M)cosfl ,^ nn . , 

£ 2 = -a 2 A- — - , (129b) 

s 2Mr ' v ' 

, sin 6* 

6 = -°* A vis- (129c) 
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For imaginary a 2 , the perturbation is no longer pure gauge. Instead it describes a change 
in the angular momentum of the black hole, given by 



j = ia 2 A. (130) 

The perturbed metric has components 

r- ■ sin 9 

h 23 = -1V2.J—, (131) 
h 22 = h 33 = 0. (132) 

We show that this is a perturbation to Kerr, written in a radiation gauge, by exhibiting 
a gauge transformation to a gauge associated with Boyer-Lindquist coordinates. Denote 
by g^"(M, J) the Kerr metric of mass M and angular momentum J. In Boyer-Lindquist 
coordinates, a perturbation h^" = -^-g^"[M, J(C)]^=o (with J(0) = 0) has as its only 
nonzero component 

^ rr = -2^sin0; (133) 
the corresponding nonzero components along the Kinnersley tetrad are 

h^ = 2h^ = -iV2j^. (134) 

With the gauge vector £ Q whose nonzero components are 



f 3 = —f 4 = - _ 
V2M 



r ( 2M 
m l 



2M V r 



(135) 



the perturbed metric h a/3 = h^l r + £^g a p is the metric of Eq. (11321) . as claimed. 

The perturbation associated with b\ is a change in the mass of the black hole of mag- 
nitude M = 3MAb\. In the radiation gauge, metric and Weyl tensor perturbations have 
components 

h 22 = -2hA, (136a) 

h 23 = 0, (136b) 

, 2(l + cos 2 #) 
^33 = M . 2a , (136c) 
sin 6 

ipo = fa = 1P2 = ^3 = ^4 = 0. (137a) 
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To verify that this is a perturbation of the Schwarzschild mass, we exhibit a gauge trans- 
formation to a Schwarzschild gauge, a gauge in which h u ^ ' \ 



h --± 



1 2M(Q 



\ r J r 

2M r 4 

-, with all other coordinate components vanish- 



r A 

ing. The corresponding nonzero tetrad components of h a p in this gauge are 

4M r 4 , M , s 

h u = — , h 22 = . 138 

r A z r 

The gauge transformation that yields a perturbed metric of this form is generated by the 
vector £ Q with components 

& = A £ [r 2 - 4Mr - 8M 2 - t(r - 2M) + 8M(r - 2M) ln(r - 2M)] , (139a) 

£ 2 = b x A [3r 2 - 4Mr - 8M 2 - i(r - 2M) + 8M(r - 2M) ln(r - 2M)] , (139b) 
£ 3 = -b ± AV2r cot9. (139c) 

The corresponding gauge perturbation is given by 

12Mr 3 

(%?) n = hA—^, (140a) 
(%?) 22 = -M — , (140b) 

r 

= M 2(l+ r/ g) - (woo 

sin v 

Finally, the gauge-transformed metric obtained by subtracting (11401) from (11361) has the form 
(fTHgjl . with M = 3MAb x , as claimed. 

The perturbation associated with b 2 is again a change in the mass of the black hole, in 
this case with M = Ab 2 . In the radiation gauge, it has the form 

h 22 = -M 2(r ~ M) , (141a) 

. , V2(r-2M)cos6 

his = -b 2 A — V 2 ' , (141b) 

r z sin 6 

, . 2(1 + cos 2 9) 

h 33 = b 2 A^ . 2Q , (141c) 

r sin 6 



00 

V>i 

"02 



03 = 04 = 0, 

Sv^Mcos^ 



b 2 A 



r 4 sin 6 1 
r - 3M 



(142a) 
(142b) 

(142c) 
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A gauge transformation to the Schwarzschild gauge (11381) is in this case generated by the 
gauge vector 

V - 2M S " 



6 

e 3 



A 



A 



2r 2 

b 2 AV2cot6 



r + 2M - 2(r - 2M) In 
,r-2M 

l + 21n 



sin^ 



sin 9 



(143a) 

(143b) 
(143c) 



for which 



(£?fiOii 

(£^)22 
(£^)33 



6 2 A 



A 5 "' 



& 2 A 



v/2(r - 2M) cos# 
r 2 sin 6 1 
2(l + cos 2 fl) 

r sin 2 



(144a) 
(144b) 

(144c) 
(144d) 



Subtracting fTll4|) from ([131)1 yields the form ffT38l) . with M = A6 2 . 

The remaining perturbations are all pure gauge. For each parameter, we list below the 
components of the associated metrics, the components of the perturbed Weyl tensor, and 
the nonzero components of the gauge vector. 



h 2 2 
h 3 3 



h 23 = 0, 
4cos6> 

c x A 



sin'fl 



2 a ' 



V> = Ipl = V>2 = ^3 = ^4 = 0, 

\/2> 



£3 = -cxA 



sin 9 



c 2 : 



^22 
^23 
^33 



V2(r- 2M) 



c 2 A 



r 2 sin 9 
4 cos 



r sin 2 ' 



(145a) 
(145b) 

(146a) 
(147) 

(148a) 
(148b) 
(148c) 
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IpO = 1p2 = 1p3 = 1p4 = 0, 

3V2M 



1> 



-c 2 A 



r 4 sin 6 ' 



(149a) 
(149b) 



fi = c 2 A 2 In ^cot Q , (150a) 

6 = c 2 A ^ In ^cot , (150b) 

6 = -^^1. (150c) 
sine* 



dr. 



h 22 = h 23 = 0, (151a) 
12M 



h 33 = -d^A—^ (151b) 



= i>i = ^2 = ^3 = ^4 = 0, (152a) 
& = diA3V2Mr(cot0 + i0sin0). (153) 



d 2 : 



h 22 = h 23 = h 33 = 0, (154) 
= i>\ = i>2 = ^3 = ^4 = 0. (155) 

Fields corresponding to ^/p 

The perturbed metric associated with \l/p has the nonzero components 

h 22 = -(5 + 2P)(5 + 4P)^p-(5 + 2p)(5 + 4p)^p 

= A— \(v 2 + tl)M 2 + t 2 x 2 — M 4 — 4M 2 tt cos + 2M 4 cos 2 0l , (156a) 
r 2 R 

h 23 = (D + p)(5 + Af3)^ P + (D + p)(5 + 4P)V P 

^/2 

= A— -f— - { (2t - M)rt M - (r - M)t 3 + (r - 2M) [rM 2 + (r - M)r 2 ] cos 
r it sin 

-2(2t - M)r M 2 cos 2 + 2M 4 cos 3 0} , (156b) 

fr 33 = p-p)(L> + 3p)*p + (L>-p)(L> + 3p)*p 
2 

= -A =- {tt 2 , - M 2 r - M(2r 2 + t 2 ,) + M 3 + 2r (M 2 + 2Mt - r 2 - t 2 ,) cosfl 

rR sm z 9 l 

[(3t - M)r 2 + (r - M)M 2 ] cos 2 - 2r M 2 cos 3 0} . (156c) 
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We now compute the change in the black-hole area for the perturbed metric associated 
with \I/p. Because the perturbed solution is static, the horizon is again a Killing horizon. 
Because the identification of perturbed and unperturbed spacetimes has been chosen to keep 
t a as the Killing vector of the perturbed spacetime, the horizon is outermost set of points 
where t a is null, where t a t a = g u = 0. Then, using h u = h 22 (this is a consequence of the 
gauge, h a pl@ = 0, and the definition of the tetrad) , we obtain the perturbation r in the 
horizon radius at fixed 9 and (f> by writing 



2M 



2M 

h 22\r=2M - ~ 



r + h 2 



22| r= 2M 
r=2M 



= -^f + A[r -M(l+cos9)}, (157) 

whence 

r = -2 AM[r - M(l + cos 9)}. (158) 

The change in the area of the horizon has, in general, contributions from the change r in 
its position and the change in the element of area, — \/^g d9d<p = \f*g (h§§ + h^i)d9d(j). In 

CIC, 

the radiation gauge, however, /i 34 = \{h^ + h^i) = 0, and the change in the horizon area is 
given by 

A-horizon = J r horizorA^ = ^ J ^ 'horizon? 'd£l (159) 

= -32nAM 2 (r - M). (160) 

The local first law of black hole thermodynamcs relates the change in the Komar mass, 
Mk = — (87r) _1 / W^t^dSap, of a Killing horizon to the change in its area: 
1 

M = 777^ — rjAhorizon- By Eq. (11601) . the particular solution associated with \l/p changes the 
mass of the horizon by 

M = -AM(r - M). (161) 

We can thus keep area and mass of the horizon constant by adding to h a p[^Ip} the / = 
metric perturbation H138[) associated with a change 

M = AM(r - M) (162) 

of the Schwarzschild mass. 

Note that instead of computing the change in area, we could have directly computed the 
change in the horizon's Komar mass. 
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As mentioned earlier, in addition to the singular radiation-gauge forms of the mass pertur- 
bation associated with b\ and b 2 , there is a radiation gauge in which the mass perturbation 

is smooth on the exterior Schwarzschild spacetime. Choosing a gauge vector £ a with com- 

2M A / r \ r 2 

ponents £ 2 = ~ hi ( — — — 1 ) , £1 = 2— £ 2 , we obtain from Eq. (11621) a perturbed 

r r 2 \2M /A 

«? 2 
metric h a p = fy^g — ^S,9a/3 whose only nonzero component is h 2 2 = — • 



F. Singularity 

For the static charge of the last section and for a static mass on a flat background (treated 
in the next section) the perturbed electromagnetic and gravitational fields are singular only 
at the position of the particle, although any choice of radiation gauge yields components 
of A a and h a p, respectively, that are singular on (at least) a ray whose endpoint is the 
particle. For a static mass in a Schwarzschild background, we have seen that the perturbed 
gauge-invariant components of the Weyl tensor are also nonsingular except at the particle. 
This is surprising, because the metric perturbation must have an additional singularity that 
can be regarded as a string or strut supporting the particle. An initially static mass with 
no external support will fall inward: If the linearized Einstein equations are satisfied for a 
source of the form T a/3 = pu a u^ \ with u a along the timelike Killing vector that Lie derives 
Sp, then the perturbed Bianchi identities imply 

= V p f a P = pu p V p u a . (163) 

It follows that the particle moves along a geodesic of the background geometry, u^Vpu a = 0, 
contradicting the assumption that the particle is static. We will find that the perturbed 
geometry has a conical singularity, and it can be chosen to run along a radial line from the 
particle to infinity. There is also, associated with any radiation gauge, a singularity of the 
perturbed metric that, like that of the vector potential in the last section, lies along a ray 
with endpoint at the particle and that can be removed by a gauge transformation. 

We now show that one can choose the parameters ai, . . . 6,2 to make the components of 
the metric finite either on a radial ray above the particle (r > r ) or below the particle 
(r < r ), but not both. If one departs from a radiation gauge, one can make the metric 
finite above and below the particle, but one cannot avoid a conical singularity. 

Although the most efficient way to obtain necessary and sufficient conditions for smooth- 
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ness of the perturbed metric is again to use the Hertz potential, in the case at hand the 
perturbed metric will not be smooth. We can at best find a gauge in which its components 
are finite, and, for clarity of presentation, we will use the direct expansion of the components 
of h a p. 

The components of the perturbed metric near 9 = 0, are given for r < r by 
4 

h 33 = A— [-2M 2 a t + b 2 + c 2 + r 2 + (61 + d - 3Md\ - r )r] 
2 

—A— {-2MV + b 2 + c 2 + r 2 

+ [3M(-ai + 3di + 1) + 61 + ci - 3M^! - r ]r} + 0(# 2 ), (164a) 
h 23 = ^fr-2M) (2M 2 0l - b 2 - c 2 - rl) + 0(61°); (164b) 

and for r > r by 

4 

^33 = [-2M 2 ai + b 2 + c 2 -r 2 + (b x + a - 3Md 1 + r )r] 

2 

-A— {-2M 2 ai + 6 2 + c 2 - r 2 
3r 1 

+ [3M(-ai + 3di - 1) + 61 + ci - 3Mdi + r ]r} + 0(# 2 ), (165a) 
/, 23 = A ^(r- 2M) ^ _ &2 _ C2 + r 2) + 0{e *y (165b) 

Then h a p is smooth at = if and only if the coefficients of 9~ 2 and 9° in /i 33 and the 
coefficient of 9~ x in h 23 vanish. The corresponding conditions on the parameters are 

2M 2 a x — b 2 — c 2 = Tq, b\ + C\ — 3Md\ = r , a\ — 3d\ = 1, r < r ; (166a) 
2M 2 ai - b 2 - c 2 = -rl, 61 + ci - 3Mcii = -r , ai - 3di = -1, r > r . (166b) 

Because the conditions for r < r are not consistent with those for r > r , one cannot find 
an single radiation gauge smooth at 9 = both above and below the particle. 
Near 9 = ir, the perturbed metric has components 

h 33 = A — 4 — [2MV + 6 2 - c 2 - (r - 2M) 2 + (61 - ci - 3M^ - r + 2M)r] 
r(7r — yj^ 

2 

-A— {2M 2 fll + b 2 - c 2 - (r - 2M) 2 + [61 -c 1 - 3Md l - r + 2M + 3M(ai + 3di - l)]r} 
+0 [(tt - #) 2 ] , (167a) 

^ 3 = A ^^y [ 2M ' ai + & 2 - c 2 - (r - 2M) 2 ] + 0[(tt - 0)]. (167b) 
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Then h a p is smooth at 9 = n if and only if the coefficients of (it — 6)~ 2 and (tt — 9)° 
in /i33 and the coefficient of (tt — O)^ 1 in h 2 3 vanish. The corresponding conditions on the 
parameters are 



2MV + b 2 - c 2 = (r - 2M) 2 , b x - a - 3Md 1 = (r - 2M), 
ai + 3c?i = 1. 



(168) 



Requiring the perturbed metric to be either 

(I) smooth everywhere except 9 = 0, r > r , or 

(II) smooth everywhere except 8 = 0, r < tq 

uniquely determines the subset of parameters {a 1; b x , b 2 , c±, c 2 , di}. For the two cases, these 
parameters - the unique solutions to Eqs. (11681) and to either Eqs. (I166ap or (1166bj) - are 
listed in Table [TV] The parameter a 2 does not affect the singular structure of h a p and is 



TABLE IV: Singularity location for two choices of parameters 



Parameter 


(I) Singularity at 


(II) Singularity at 




= 0, r > r 


= 0, r < r 


ai 


1 













h 


(r - M) 





b 2 


-2M(r - M) 


rl - 2Mr + 2M 2 


Cl 


M 


-r + M 




-r (r - 2M) 


2M(r - M) 


di 





i 

3 



arbitrary; for definiteness it is set to zero in the table. The trivial parameter d 2 does not 
alter the metric and is not listed. 

The perturbed metric h^l corresponding to the parameters of the first column is, in 
particular, smooth for r < r$. The nonzero values of the parameters b\ and b 2 correspond 
to a perturbation of the Schwarzschild mass given by 

M = 3AMb 1 + Ab 2 = AM(r - M). (169) 
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This is exactly the change in mass (11621) needed to cancel the change in black hole area and 
mass of the particular solution. Smoothness of a radiation gauge in the region 2M < r < tq 
between the horizon and the particle thus requires that the perturbation involve no change 
in the black hole mass. Just as the vector potential in a radiation gauge is smooth on a 
sphere only if the sphere encloses no perturbed charge, our result suggests that the perturbed 
metric can be smooth on a sphere only if the sphere encloses no perturbed mass. 

We thus expect that for the perturbed metric h^P , smooth for r > r (and asymptoti- 
cally flat), the asymptotic mass must vanish, and direct calculation verifies this. For a flat 
background, we will see in the next section that one can obtain from h^P a perturbed met- 
ric that agrees up to gauge with h^p by adding an I = part, a perturbed Schwarzschild 
solution centered at the origin that restores the total asymptotic mass. For a Schwarzschild 
background, however, there is a second physical difference between solutions (I) and (II): 
The difference in the values of ai means that hrP and h^P differ by a perturbed C metric. 



The solution that agrees up to gauge with h ( P has the form h l *P + h s a f warzschad + h^ metric 



*/9 allu lb af3 

*/3 has the form n aj3 -r u ap -r 
The C- metric term has a conical singularity on the axis; adding it in the gauge of 26j], and 
adding the Schwarzschild perturbation in a Schwarzschild gauge yields a perturbed metric 
that is finite but not smooth on the axis. 



V. FLAT SPACE 

When the mass M of the background Schwarzschild geometry is set to zero, the per- 
turbation is just the gravitational field of a point mass linearized about a flat background: 
the linearized Schwarzschild solution. Like the electromagnetic field of a point charge in a 
Schwarzschild background, the perturbed geometry is now singular only at the position of 
the particle, and any other singularity in the tensor h a p is an artifact of the choice of gauge. 

The origin of coordinates inherited from the M/0 case is displaced from the particle by a 

distance r ; and the inherited tetrad is aligned with inward and outward radial null geodesies 

associated with the spatial origin, not with the particle. With respect to a tetrad associated 

with the null geodesies that start or end at the particle, the perturbation is algebraically 

special, having as the only nonzero component of its perturbed Weyl tensor ip2 = with 

R 6 

R the distance to the particle. With respect to our null tetrad, however, the perturbation 
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is not algebraically special: Its components are 

m 



V>o = -3-^r 2 sin 2 £, (170a) 
3 tn 

ih. = -^=-^r (r -r cos9)wi9, (170b) 

^ = -^(^-l^ej, (i70c) 

^3 = -^^-^r o (r-r o cos0)siii0, (170d) 

^4 = -^r 2 sin 2 0, (170e) 



where 



R = Jr 2 + r 2 - 2rr cos9. (171) 



(This coordinate expression for the distance R to the particle is the M = limit of the 
length R of the last section.) In contrast to perturbations of Kerr and Schwarzschild, where 
only ip and ip^ are gauge invariant, the full Weyl tensor is gauge invariant, because the Weyl 
tensor of the flat background vanishes. The corresponding components along the null tetrad 
l a ,n a , m' a differ from these only by a phase and are manifestly smooth except at R — 0, the 
position of the particle, and at r = 0, where the tetrad is singular. 

^ = - 3 |ro sin2 ^ = -^i* + n/) 2 > ( 172a ) 

^ = 7fS r ° (r_r V )sin ^ = ^f r f(r- r fz)(x + l y), (172b) 



R 2 - 3 {f 2 (x 2 + y 2 ) , (172c) 



^ = -^^r (r-r cose) 8 m^=- i ^^(r-^)(x-iy), (172d) 

K = -\^™ 2 = -l^{*-iv) 2 - (172e) 

The components ^[ are smooth at 6 = n, although the tetrad is not; this is due to the fact 
that the tetrad l a ,n a ,e~' l ^m a is smooth at 9 = ir, and the components of the Weyl tensor 
associated with that tetrad are simply 
The perturbed mass is given by 

lim (— r s ip2) = lim r 3 ^- ( R 2 — -rQsin 2 # J 

r^oo ' r^oo it V 2 / 

= m. (173) 
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We will again find the perturbed metric h a p in a radiation gauge, choosing, within the 
family of radiation gauges, one for which the metric components are finite except on the 
9 = axis for r > r or for r < r . As was the case for the vector potential of the smooth 
electromagnetic field of a point particle, one cannot globally choose parameters a 1; . . . d 2 for 
which the metric is smooth and regular everywhere outside the particle. 

As in the previous sections, if the Hertz potential associated with the smooth tetrad 
l a , n a , m' a is smooth at 9 = 0, then the perturbed metric associated with this tetrad is 
smooth at 9 = 0. 

From Eqs. (11221) and (11231) . a particular solution \l/p associated with the Kinnersley tetrad 
has the form 



2 A R 3 
—A — -„ 

3 sin' 9 

and the general homogeneous solution has the form 



*p = -0^—2^ ( 174 ) 



V H = A — \( a i cos3 8 - 3q i cos # + 6di)r 3 + 3(&i cos 2 9 + hi + 2c x cos 9)r 2 
3 sin 9 

+ (a 2 cos 3 9 - 3a 2 cos6» + Qd 2 ) + 3(b 2 cos 2 9 + b 2 + 2c 2 cos 9)r] . (175) 

The Hertz potentials associated with tetrads smooth at 9 = and 9 = tt are then ^e ±2l< ^ = 
(V P + ^ H )e ±2i<t> . 

We will show that one can make the corresponding metric h a p smooth at 9 = tt, by a 
choice of parameters that correspond to gauge transformations. One can simultaneously 
choose parameters to make h a p smooth at 9 = for either r > or for r < 0. But no 
consistent choice of parameters allows the perturbed metric to be smooth both above and 
below the particle on the ray through the particle. Moreover, smoothness at 9 = above 
the particle requires a choice of parameters that alters the asymptotic mass; to recover for 
r > To the metric of a static particle, one must add to the radiation-gauge solution an / = 
part, a linearized Schwarzschild solution centered at the origin. 

We begin by describing changes in the perturbed metrics associated with the parameters 
oh , . . . , d 2 that arise for a flat background. For Kerr geometries with nonzero mass, Wald 



321 ] shows that any perturbation with ipo or ^4 = is a type D perturbation or is pure 
gauge. Because the proof relies on the fact that a gauge transformation can make ipi an d 
ip3 vanish, it fails for a flat background. There are, in fact, vacuum perturbations of flat 
space for which the gauge invariant components if>i, ip 2 and ^3 have nonzero values, although 
ipo = ipA = 0. In particular, we will see that the perturbation associated with a 2 is no longer 
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pure gauge but is a perturbation of this kind. (This can happen because the gauge vector 
ffT29D has no M = limit). 

Conversely, we will see that the perturbations associated with a\ and b 2 are now pure 
gauge, although they are C metric and mass perturbations when M/fl. The components 
of the perturbed metric and Weyl tensor associated with each of the parameters can be read 
off from the corresponding expressions in the last section by setting M to zero. Because 
each parameter corresponds to a vacuum perturbation and the background Riemann tensor 
vanishes, the perturbation is pure gauge if and only if the components of the perturbed Weyl 
tensor vanish. That is, a perturbation of flat space is pure gauge iff the perturbed Riemann 
tensor R a ^s vanishes; and, for a perturbation of flat space, we have 

-^Vrtf = <^ Rafiyd = ^ R a/ 3 = and C a p r s = 0. 

Eqs. (I126p imply that the perturbed Weyl tensor associated with a\ vanishes, and thus 
that di corresponds to a gauge perturbation, h a p = £^ a/ g. 

Writing A = A|m=o = —■, we obtain 

h 22 = —aiA2rcos6, 

h 23 = —aiA —=rsin9, 
v2 

h 33 = 0, 



(176a) 
(176b) 
(176c) 



£1 = axAhyt -r) 2 cos#, (177a) 
£ 2 = - ai A^(t -r)(£ + 3r)cosfl, (177b) 
6 = - ai A^=(t 2 -r 2 )sm6. (177c) 



Because the components (11281) of the perturbed Weyl tensor associated with a 2 are 
nonzero, the perturbation associated with a 2 is no longer pure gauge. 

The perturbation associated with bi was a change of mass proportional to M; because 
Eqs. ( 11371) imply ipi = 0, it is pure gauge when M = 0: 

h 22 = -2bxA, (178a) 

h 23 = 0, (178b) 

, ~ 2(l + cos 2 fl) 
h 33 = hA . ; , 178c 
sin 
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& = M (r-f), (179a) 
6 = M^(3r-*), (179b) 
£ 3 = -b ± A V2rcot9. (179c) 

The metric associated with b 2 remains a mass perturbation, a linearized Schwarzschild 
solution, centered at the origin with mass — mb 2 /rQ . Finally, the perturbations associated 
with Ci and c 2 are again pure gauge; and the perturbations associated with both d\ and 
d 2 now vanish (h a p = 0). To summarize: In flat space the perturbations associated with 
ai, bi, ci, c 2 are pure gauge, and those associated with d\ and d 2 vanish. 

Then one can find a smooth perturbed metric of a static particle in flat space from a 
radiation-gauge Hertz potential only if one can choose the six parameters a±, &i, Ci, c 2 , d±, d 2 
to make the perturbed metric smooth. Now h a p is smooth on the axis if and only if its 
components in a smooth basis are smooth. Because h 2 ry and hyy have angular dependence 
e imii> and e 2im ^ , they are smooth at 9 = only if h 22 = 0(6°), h 23 , = 0(6), h wy = 0(9 2 ); and 
this is true if and only if 

h 22 = 0(6°), h 23 = 0(6), h 33 = 0(6 2 ). (180) 
Similarly, the metric can be smooth at = tt only if 

h 22 = 0((9-n)°), h 23 = 0((9-n)), h 33 = 0{{6 - tt) 2 ). (181) 

These requirements impose more algebraic conditions on the parameters than can be 
simultaneously satisfied. One can, however, choose parameters for which the metric is 
smooth everywhere except on a ray that does not intersect the particle (and at the particle 
itself): That is one can find a perturbed metric from a Hertz potential ^' that is smooth 
everywhere except 9 = ir. 

As we saw with a Schwarzschild background, one can also choose parameters that make 
the perturbation regular in a radiation gauge either inside or outside the particle. Inside 
the particle, the choice is simply a gauge transformation. Any choice of parameters that 
make the perturbation regular outside the particle require zero net mass for r > tq. One 
can obtain a regular solution in a radiation gauge for the / ^ part of the perturbation; to 
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complete the solution one must add an / = part - corresponding to a mass at the origin - 
in a different gauge. 

One can find the allowed parameters either by directly expanding the metric for arbitrary 
parameter values near = and = n or, more simply, by examining the Hertz potential 
near = and = it. Now \J r ' is smooth near = (0 = tt) if and only if \1/ = O(0 2 ) 

= 0((9 - tt) 2 )). From Eqs. fll74H75j) . the parts ^ P and ^ H of the Hertz potential have 
near 9 = the forms 



-A-\r — r n | 
3 1 1 



(r - r ) 2 9 2 + -(r - r ) 2 + ^r r 



+ O(0 2 ), 



tt* = 4{ + 3rf i) r3 + 3 ( 6 i + c i) r2 + 3 ( & 2 + c 2 )r - a 2 + 3d 2 ] 0" 



+ 



-(-ai + 3rfi)r 3 - -(6i + ci)r 2 - 1(& 2 + c 2 )r + -(-a 2 + 3d 2 ) 



(182a) 



O(0 2 ). 
(182b) 



Then, for *f> = typ + we have 



* = i| I [(-1 - a x + 3rfi)r 3 + 3(r + h + Cl )r 2 + 3(-r 2 + 6 2 + c 2 )r + r 3 - a 2 + 3rf 2 ] 0" 



+ 



1 111 

-(-1 - ai + 3di)r 3 - -(r + 6 X + c x )r 2 - -(-r 2 + 6 2 + c 2 )r + -(r 3 - a 2 + 3d 2 ) 



+O(0 2 ), r>r , (183a) 
tt = 4|[(l-a 1 + 3di)r 3 + 3(-r + &i + c x )r 2 + 3(r 2 + 6 2 + c 2 )r - r 3 - a 2 + 3rf 2 ] 0~ 2 

-(1 - ai + 3rfi)r 3 - -(-r + &i + c x )r 2 - l(r 2 + 6 2 + c 2 )r + l(-r 3 - a 2 + 3rf 2 ) 



+O(0 2 ), r < r„. 



(183b) 



Similarly, for near 7r, we have 



* = i~ | [(-1 + ai + 3rfi)r 3 + 3(-r + h - Cl )r 2 + 3(-r 2 + 6 2 - c 2 )r - r 3 + a 2 + 3d 2 ] (0 - tt) 



+ a x + 3d 1 )r 3 - i(-r + 6i - c x )r 2 - l(-r 2 + 6 2 - c 2 )r + i(-r 3 + a 2 + 3d 2 ) 



+ 

6 

+0[(9-n) 2 }. 



(184) 



We thus obtain as necessary and sufficient conditions for smoothness of tt at = 0, 



a% - 3di = —I, a 2 - 3d 2 = r%, b\ + c\ = — r , b 2 + c 2 = r%, r > r ; (185a) 
ai - 3d 1 = 1, a 2 - 3d 2 = —r%, h + c x = r , 6 2 + c 2 = -r 2 ,, r < r ; (185b) 
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and at 9 = tx 



ai + 3di = 1, a 2 + 3d 2 = r%, h — ct — r , b 2 - c 2 = r%. (186) 

Because the perturbations corresponding to d% and d 2 vanish, one can satisfy the first 
two conditions in each line without changing h a p. The necessary and sufficient conditions 
for smoothness of h a p at 9 = (as can be verified directly from its components) are then 

h + ci = -r , b 2 + c 2 = rl, r > r ; (187a) 

bi + ci = r , 6 2 + c 2 = -rjj, r < r ; (187b) 

and at # = 7r, 

6i - ci = r , b 2 -c 2 = Tq. (188) 

For the two cases, these parameters - the unique solutions to Eqs. (11881) and to either 
Eqs. f TTSTaD or ffl87b| - are listed in Table ED 

As in the last section, Eqs. (11871) immediately imply that one cannot find a radiation 
gauge that can be everywhere locally obtained from a Hertz potential, for which h a p is 
simultaneously smooth at 9 = outside and inside the particle. One can choose parameters 
for which h a p smooth everywhere except 

(I) along the part of the 9 = ray below the particle (with r < ro), or 

(II) along the part of the 9 = ray above the particle (with r > ro). For the parameter choice 
that makes h a p smooth outside r = ro, however, b 2 has the nonzero value — changing 
the asymptotic mass to zero by subtracting a mass m Schwarzschild solution centered at 
the origin. As in the Schwarzschild perturbation of the last section, one can then obtain a 
smooth metric outside r = r with asymptotic mass m by adding back an / = perturbation 
in a smooth gauge - adding, for example the linearized Schwarzschild solution of mass m 
centered at the origin. 

The perturbed metric smooth everywhere except 9 = 0, r > ro is given by 

4? = "2m (J- - , (189a) 

hg ] = -V2m-^— (^-- - 4 cos ^ > ( 189b ) 
16 sin \rR r R I 



h$ = -2m- 



sin 2 fl 



i I +2 (I_^_^_) cos9 + (i_I|,w« 

K r \ r r tt r oR J \R r o 



. (189c) 
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TABLE V: Singularity location for two choices of parameters 



Parameter 


(I) Singularity at 


(II) Singularity at 




9 = 0; r > ro 


9 = 0; r < ro 


ai 
















h 







b 2 





r o 


Cl 







C2 


r o 






The perturbed metric smooth everywhere except 6 = 0, r < r is given by 



h. 



(ii) 

22 



,4P = -V2 



= -2m- 



sin 2 



1 



11 

«-; +2 '-- 



—77 ) cos0 + ( - ^ ) cos 2 f/ 



(190a) 
(190b) 
(190c) 



r r i? ri?, 

The perturbed metrics h^p and + /i^ hwarzschlld (m) must agree, up to a gauge trans- 
formation, with the perturbation one acquires from the Schwarzschild metric by translating 
the origin and linearizing about flat space. This is the case, and the explicit gauge vector 
£ a for the metric h T a/3 has components 



6 = 

6 = 
£3 = 



— m 



t — r r — r n cos 9 

1- 

r R 

m / t — 3r r — r cos 



2 In 



r — r n cos 9 + R 



sin 



2 



+ 



+ 2 In 



r cos 9 + R 



sin 



m 



2^2 sin ir 



4r „ 2LR 2 + (r - r o cos0) 2 , 
— cos H — — — 4 



r R 



(191a) 

(191b) 
(191c) 
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The resultant metric is given by 

2m 



h n = ^(2i? 2 -r 2 sin 2 #), (192a) 

h 2 2 = ^(2i? 2 -r 2 sin 2 fl), (192b) 

^33 = (192c) 

hl2 = T^l t ( 192d) 

\/2mr (r -r cos 9) sin0 

^13 = ™ , (192e) 



mr (r - r cos6») sin6» 

ft23 " ' ( ' 



VI. FINDING THE SELF-FORCE IN A RADIATION GAUGE. 



An algorithm for finding the self force on a mass moving in a background spacetime, 



the MiSaTaQuWa method, was given several years ago lOl , In this method, as noted 



21], a particle follows a geodesic of a 



by Quinn and Wald and by Detweiler and Whiting 
renormalized metric, h ren , given by 

h ien = h vct - h sing , (193) 

where h ret is the retarded field of the particle in a Lorenz gauge and h smg is its locally 
defined singular part, chosen to cancel the singular part of h mt and to give no contribution 
to the self-force. In the Quinn- Wald characterization of the method h sms is obtained in a 
neighborhood of a point P of the particle's trajectory from the field of a particle moving in 
flat space by using the exponential map to relate the linearized field in flat space to a field 
^smg j n neighborhood G f P. One demands that the image under the exponential map of 
the flat-space particle's position, 4-velocity, and acceleration agree at P with those of the 
original particle. 

Because it involves a Lorenz gauge, the method cannot take advantage of the decoupled, 
separable wave equation - the Teukolsky equation - that governs black hole perturbations. 
One could imagine constructing the singular part of the gauge- invariant perturbed Weyl 
tensor component \l/o (or ^4) from h smg and then obtaining the renormalized metric in a 



47 



radiation gauge from \1/Q 6n (or W^ n ). The radiation-gauge prescription, however, requires the 
perturbed Weyl scalar to be a solution to the linearized Einstein equation. In the forms in 
which it is given by Mino et al.flo| and by Quinn and Wald 11], h smg is not a solution to the 



Einstein equations linearized about the background Kerr spacetime; and it is not obvious 
that one can obtain a prescription for a renormalized Weyl scalar \1/Q en (or \l/4 Cn ) that satisfies 
the Teukolsky equation. 



A more recent version of the MiSaTaQuWa method, due to Detweiler and Whiting [21], 
overcomes this difficulty In their prescription, the singular part of the field h smg is redefined 
in a way that makes it a solution to the linearized Einstein equations, defined in a normal 
neighborhood of the particle. As a result, one can compute a gauge-invariant renormalized 
Weyl scalar in the following way: 

The perturbed Riemann tensor is given in terms of h a @ by 

From Eq. ([7j), i/jq has the form, 

^ = -l a m p Pm s C aPlS = -FmPrnfRafrs = O aP h aP , (195) 
where O a/3 is the operator 

af3 = l(i*i0 m i m s + m a m Ppi s - l a m 0p m s - / Q m /3 m 7 / <5 )V 7 V 5 . (196) 

In particular, from the singular part of the perturbed metric, h s ^ s , one can define the singular 
part ^o mg °f the perturbed Weyl tensor by 

Then ipQ en is given by 

VC = Vo et ~ rltf*, (198) 

with 

Because h iet and h smg have the same source, so do ip Tet and ip smg . Their difference, ip ren is 
thus a sourcefree solution to the Teukolsky equation. That implies that one can construct a 
nonsingular h ren in a radiation gauge from \|> ren . 
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In a Lorenz gauge, we have 



where h a/3 = h a/3 q a/3 h and where the Green's function has the form 

2 

sin 1 / I 

G aW5>( x i x> ) = yU a py S ,(x, x)8(a) + -V a ^> s >(x, X )6(o"). (201) 



dere U and V are nonsingular bitensors, whose form is given, for example, in [2l[ (see also 
l|). The step function 0(er) has value 1 when a > - when x and x' are spacelike separated; 

and it vanishes for a < 0. Thus G^/jytf' ( x > x ') nas peculiar causal behavior: It is nonzero 

only when x and x' are not timelike related. 
A Green's function for \I/Q ing is then given by 

Gff{x, x>) = (<y# - \g^O\} G s ^, s ,(x, x'). (202) 

Having computed G s ^(x,x'), one can obtain ipQ ng from the expression 

^ ing (;r) = [ G^s'&x'W^'ix'^V. (203) 



The explicit forms of U and V have been computed in terms of the geodesic distance from 
the particle, to the order needed to find the self-force (although the order needed for the 
prescription given here may be different). 

From vpQ en one recovers h™a, using the CCK equations (I102I) or H106[) in a Kerr background 
to find the potential ^ and using Eq. (|93|) to find tf en in terms of Eq. fl 102j) can be 
solved by four radial integrations for each angular harmonic if one works in the frequency 
domain. The operator D, restricted to an m, ui subspace corresponding to and t dependence 
e i ( m *- w *), has the form 

\ma — (r 2 + a 2 )cu] 
A™ = d r + i [ - K — (204) 

Then \l/ ren is found from 

^ = D L^Z- (205) 
Because h™p is a solution to the linearized vacuum Einstein equations, it is determined by 
■00 (or -04) up to an algebraically special vacuum perturbation. 

The method of finding the self-force in a radiation gauge may now be summarized as 
follows: 
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1. Compute ^Q ing either from the known form of h sms as outlined above, or, with addi- 
tional insight, directly from the Teukolsky equation. 

2. Express ipQ ag in terms of spin- weighted spheroidal harmonics. 

3. Compute -0Q et from the Teukolsky equation. 

4. Write ipo en — V^ 6 * — regularizing by a cutoff in angular harmonics, and finding 
the limit, to desired accuracy, as the cutoff harmonic is increased. 

5. Find the potential ^ either from ip r ^ n , using Eq. (11061) . or from ijjQ en by 4 radial inte- 
grations of Eq. (12051) . 

6. Find the perturbed metric h a p in a radiation gauge by taking derivatives of This 
leaves spin-weight and 1 parts of the perturbation undetermined. 

7. Obtain the spin-weight and 1 parts of h a p by fixing the area and angular velocity 
of the perturbed black hole and using jump conditions across the particle of the spin- 
weight and 1 parts of the perturbed Einstein equation. 5 

8. Compute the self-force from the perturbed geodesic equation. For the perturbed metric 
9a/3 + (h aj 3, with V the covariant derivative of g a p and u a normalized by g a p, the 
geodesic equation has, to 0((), the form 

u ■ Vu a = -C(g a5 + u a u s ) (v p }q? - \v s hf^ = f . (206) 

With the self-force f a computed, one must find the particle trajectory as a self-consistent 
solution to Eq. (12061) . The emitted radiation can be found from ip^. 

Because the perturbed Weyl tensor involves two derivatives of the perturbed metric, it is 
possible that the renormalization program outlined here will require a Hadamard expansion 
that is two orders higher in the separation between source and field point than the expansion 
used in a Lorenz gauge. We suspect, however, that the extra orders can be avoided by using 
Eq. (12051) to find the Hertz potential. The argument is that each radial integration changes 
by one factor of R the highest power of 1/R and thereby reduces the order of the singularity. 

5 Using the jump conditions of the spin-weight and 1 parts of the Einstein equation was suggested to us 
by Bernard Whiting. 
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Calculations now underway in a Schwarzschild background will decide the issue and may 
show whether the method can feasibly be used to obtain generic orbits in a Kerr background. 
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APPENDIX A: SMOOTHNESS IN TERMS OF A HERTZ POTENTIAL 

In the CCK formalism one writes the vector potential A a and the perturbed metric h a p 
in terms of Hertz potentials $ and The components of A a and h a p, given by Eqs. fl53|) 
and (11001) . have in each case the form (or where £ is a sum of covariant derivatives 
along the basis vectors and Christoffel symbols (spin coefficients) of the basis. Because 
the Christoffel symbols of a smooth basis are smooth, smoothness of the Hertz potential 
$ associated with a smooth basis {e^} is a sufficient condition for smoothness of A a ; and 
smoothness of the Hertz potential ^ associated with a smooth basis is a sufficient condition 
for smoothness of h a p. 

We next relate this condition to a condition on the Hertz potential associated with the 
Kinnersley basis, by showing that the Hertz potentials $ ± and \l/ ± associated with the 
smooth bases {e^} differ only by a phase from the Hertz potentials $ and In addition, 
the metrics corresponding to the two sets of Hertz potentials are identical. 

These results are immediate from the following proposition, proved here for a 
Schwarzschild geometry, but presumably correct for type D vacuum spacetimes. 
Proposition, (i) Let ^ satisfy the sourcefree Bardeen-Press (Teukolsky) equation for 
spin- weight s = 2, namely T^l = 0. Then ^/^ = ^/e T2l< ^ satisfies the same sourcefree 
equation, = 0, and /i Q/3 [^ ± ] = h a p[^\. 

(ii) Let $ satisfy the sourcefree Teukolsky equation for spin- weight s = 1, T$ = 0. Then 
$ ± = $e T ^ satisfies the same sourcefree equation, T5> ± = 0, and A Q [ ( I )± ] = A a [<&]. 
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Proof. The proof essentially follows from the fact that expressions arising in the NP and 
CCK formalism have definite spin- weight. A Lorentz tranformation of the basis that changes 
only the phase of m a changes the phase of each expression only by a phase corresponding 
to the spin weight of that expression. The detailed verification is as follows. The operators 
and spin coefficients associated with the bases {e^} are related to the corresponding objects 
of the Kinnersley basis by 

D ± = D, A ± = A, 5 ± = e^S, 5* = e^H (Ala) 

/3± = _ a ± = p - 1 ■ - e * 7 ± = 7 , p± = p, ^ = n. (Alb) 
2v2rsm0 

For a quantity 77 having spin-weight 2, we have 77 =l= = e ±2l( ^rj, and from Eqs. flAlj) above we 
obtain the relations 

(D + np) ± 7y ± = e ±2 ^(L> + np)r), any n, (A2a) 

(A + 2 7 + pf 77+ = e ±2i * ( A + 2 7 + p) rj, (A2b) 

(5 + 4/3)^ = 6^(5 + 4/3)77, (A2c) 

Similarly, a quantity rj of spin-weight 1 satisfies the relations 

D ± 7 1 ± = e^Dr], (A3a) 
(5 + 2/3)V = (5 + 2/3)77, (A3b) 
(5 - 2 / 3) ± r ? ± = e ±2 ^(5 - 2/3)ij. (A3c) 



From the form of the Teukolsky spin-2 and spin-1 operators of Eqs. (110 II) and (1541) . 

T s=2 = ( A + 2 7 + fi) (D + 3p) - (5 - 2/3) (6 + 4/3) - 3^ 2 , ( A4) 

T s=1 = (A + jj,)(D + p) - 5(8 + 2/3), (A5) 

and Eqs. ( IA21) and ( lA3~l) . we have the claimed relations T ± \l> ± = T*, T ± $ ± = T$. 

From the form (llOOp of h a /3, Eqs. (1A2|) imply that the components /io,^^ 111 ] along the 
basis {e^} are related to the components of /i Ql g[\l/] along the Kinnersley basis by 

fc 22 [tf±] = fr 22 [*], h 23 = er^M*]. M*1 = e ±ai *M*]- (A6) 
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But these are just the components of /i a( g[^] in the basis {e^}, implying /i^t^] = ft-o^f 1 ^], 
as claimed. 

Similarly, from the form (J53D for A a , Eqs. ( [Ml) imply A ct [$ ± ] = A a [$\. □ 

APPENDIX B: <j> 2 AND ^4 IN CLOSED FORM FOR STATIC PARTICLES. 

We first show that the expression (jHJ) for 02 as a series has the sum f j4"5l) and then turn 
to the gravitational case, showing that the series (1901) has the sum ( 1911 . 

Noting that the sum (JHJ) is axisymmetric and using Eq. ffTSl) . we can write the sum in 
the form, 



. 47re(AA 

>2 ^ 



1 1/2 



2V - 2Mr , ro ^wTW^ Pl ^) Q; ^)- lM0 ' mo) < B1) 
= 5 g<<rV' (£)*<MW>,o) (B2) 

To sum the series HB3I) . we will use the relation 



CO 

S = = 5^(21 + l)fi(a)Q,(6)fi(x), b > a > 1, (B4) 

a/gt + tr — 2aox — (1 — x 2 ) " 

which follows immediately from the series expression for the Laplacian Green's function, 
1/R, written in prolate spheroidal coordinates [33J. That is, the distance R between points 



with prolate spheroidal coordinates (r,9,(f>) and (ri,#i,0i) is given by 



R 2 = r 2 +rf — 2rri cos 9 cos 9\ — 2^/r 2 — e 2 \j r\ — e 2 sin 9 sin 9\ cos(0 — <pi)— £ 2 (sin 2 #+sin 2 (9i), 

(B5) 

where r\, r > e; and the Green's function is given for ri > r by the series 



1 _ 1 

R ~ 1 



oo 

^(2/ + 1)[P, (£) (^) Pi(cos0)P,(cos0i) + 



Z=0 

2 f^Kf^f} P ™(~) ^(7) ^ m (cos^)Pr(cos^) cosm(0 - &)] 

m=l k \ / " ^ 



.(B6) 



7" f\ 

Setting 9\ = 0, - = a, — = 6, and x = cos 9 in Eqs. ( 1B5I) and ( 1B6[) . we obtain Eq. ( 1B4I) . 
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We next use the standard identities 34], 

Pi(x 



[x 



(2l + l)x 

1) ^ 

dx 



Qi(x 

P t (x 

Qi{x 



= (1+1) 
(1+1) 



Pl+lfl)' 


+ i 


Pi-i(x)' 
Qi-i(x)_ 




— X 


~Pi(x)~ 


1 


_Qi+i(x)_ 


.Qi(x). 





to show the relation 

Si = (— ab + x)S + a 



(a 2 - 1)V^ _ 1)a / 2 £ ^±^pi (fl) gi (6) ^ (s 



(B7a) 
(B7b) 

(B8a) 
(B8b) 



(B9) 



First, from the definition (IB4I) of 5" and Eq. (lB7a[) . we obtain 

oo 

= ^^_i(a)g,_i(6) + + l)iVf.i(a)Qj+i(6)]^(^)- 
Now, substituting Eq. (IB9I) in Eq. (IB8a[) . we have 

oo 

S-i = - 2^(21 + l)aP(a)&Q*(6) - /fl_ 1 (a)Q,_ 1 (6) - (Z + l)fl +1 (o)Q I+1 (6)]^(x) + a 

oo 

- ^[(2/ + l)a^(o)6Q,(6) - ZfU(a)Qi-i(&) - (I + l)P z+1 (a)Q z+1 (&)]P(x), (BIO) 



where the relations P = 1> -Pi ( a ) = a ; Qiip) = bQ (b) — 1, are used to eliminate the / = 
terms. Again from Eq. (lB7al) . we obtain 



Si 



i=i 



E ( 2/ + 1 ) 



(1+1) , s 1 

-Pi+xia) + , nl , ^ Pi-i (a) 



(21 + iy 



(21 + iy 



(i + i) 

(21 + 1) 



Ql+i(b) 



I 



(21 + 1) 



Qi-i(b) 



lPi-i(a)Qi-i(b) - (I + l)P +1 (a)Q m (fc) 



- P/-i(»)Q /+ ,(M - P / _ J ( ( oQ,- J (^) + P,_,(«)Q/-i(fc)] Pf.r) 



oo 



2/ + 1 



z=l 



2/ + 1 



(Bll) 



Using the definitions, 

PT(x) = (* 2 - i) m/2 ^(*), QT{x) = (x 2 - l ) ml2 ^M*) 



\x\ > 1, 



and the identities ( lB7bl) and ( 1B7al) . we can write 



(B12) 



(B13) 
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Using this relation and the analogous relation for Q] (b) to replace the bracketed expressions 
in Eq. flBTTl) yields Eq. flB8bl) . as claimed. 



Finally, from Eqs. (1B3P and (1B8bl) . with a = — , b = — , x = cos9, we obtain the closed 



form expression 



for 



> 2= e 't A ° )V2 § 



2 v / 2A'f 



M tt -M 2 cos# 



(AAo) 1 ^ 



i? 



A A sin 9 



(B14) 
(B15) 



2V2r r 2 R 3 

Next, we turn to the gravitational series. Again we use axisymmetry and Eq. (fISl) to 
write the sum (190]) in the form 

1/2 OO r,. n,-|1/2 



1p 4 



vrm A A 



M r 4 r 
vrm A A 



1=2 
1/2 x 

^8 2 V 



(*-2)! 
(Z + 2)! 



r 



2^0^,0)^(0.0) 



M r 4 r " ^( / + 2 ) ! 



M 



M 



m AA 



1/2 x 




4M r 4 r 



8 2 E(^i)^(S)o?(|)fl(c OS , 



In this case, we wish to show 



(B16) 



a 2 + b 2 + a 2 b 2 - 1 - Aabx + 2x 2 )S - (a 2 + 1)6 - a(a 2 - 3)x 



(a 2 - 



1)(^ 2 -1)E 



(2/ + 1) 



^ (I + 2)(l + 1)1(1-1) 
Using Eq. ( 1B9I) . Eq. ( 1B7al) . and relabeling the indices, we obtain 



if(a)Q 2 (6)P^)- 



x 



1=0 



(2/-1)" 

■ (* + !)(* + 2) 
(2/ + 3) 

Similarly, from Eq. flB7al) . we have 
(/ + !)(/ + 2) 



(Z + l) 2 I 2 



(21 + 3) (21 - 1) 



Pi(a)Qi(b) 



Pi + 2(a)Qi +2 (b)\Pi(x). 



a 2 P t (a) 



(2/ + l)(2/ + 3) 
1(1-1) 



Pi+2(a) ^ 
Pi-2 (a) 



/ 2 



(2/ + l)(2Z + 3) (2Z + 1)(2Z-1) 



(B17) 



(B18) 



P,.(a) 



(2Z + 1)(2/ - 1) 
together with the corresponding expression for b 2 Qi(b). 



(B19) 
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As in the electromagnetic case, the I = 0, 1 pieces of S2 that do not multiply S are chosen 
to cancel the / = 0, 1 pieces of the series, so that 5*2 can be rewritten as 

00 

S 2 = J^(a 2 + b 2 + a 2 b 2 - 1 - Abx + 2x 2 )(2l + l)P(a)Q ; (&)P(x) (B20) 

1=2 

00 

= Y,( A + Bx + 2x 2 )(2l + 1)^,(0)^(6)^,(2:). (B21) 
1=2 

We use Eq. (IB19I) (and the corresponding expression for b 2 Qi(b)) to write the part of 
Eq. (lB2ip involving A in the form 



J2 A (2l + l)P l (a)Q l {b)P l (x) 
1=2 

(/ + !)(/ + 2) 



E 



1=2 

+Pi{a 
1 

+ 2l + 1 

x 



(21 + 3) 

(/ + !)(/ + 2) 
(21 + 3) 
(/ + !)(/ + 2) 

(21 + 3) 
(/ + !)(/ + 2) 



Pi+*W + 

Qw(b) + 
Pi+2(a) + 



(l + l) 2 I 2 



(21 + 3) (2/-l)_ 
(/ + 1) 2 + I 2 



P(a) 



*(*-!) 
(21 - 1) 



Pi- 2 (a)}Qi{b) 



(21 + 3) (21 - 1) 
' (I + I) 5 



+ 



I 2 



+ 



(2/ + l)(2/ + 3) 



Qi + 2(b) + 



(21 + 3) (21 - 1) 
(^ + 1) 2 



P(a) + S^4rP,_ 2 (o) 



(2/-1) 



+ 



_(2/ + l)(2/ + 3) (2/+ 1)(2Z — 1)_ 
(2l+l)P l (a)Q l (b) )Pi(x). 



Qi(b) 



(21 + l)(2l - 1) 

We next use Eq. flBTal) to rewrite Ya=2 Bx ( 21 + l)P,(a)Q,(6)P,(a;) as 



(B22) 



B x(2l+l)P l (a)Q l (b)P l ( 
I 



x 



1=2 



00 

Z=2 



(2/ - 1) 
(J + 2) 



P,(a) + ^-i-P,_ 2 (o) 



(21 - 1) 
(i + 2) 



Qi(&)+ttHtGi-2(&) 



Qw(b) + 



2/ - 1 
l + l 



Qi(b) 



Pz(x).(B23) 



(2/ + 3) ^ v ' 2/ + 3 v 'JL(2/ + 3)^'^ V ^ ' 2/ + 3 
Substituting in Eq. flB2Tj) the expressions from Eqs. (IB22I) . (1B23I) . and (1B18I) . we find 

(l + 2)(l + 1)1(1-1) 



S-2 



(71 - 



1=2 



(21- 1) 2 (2/ + 1)(2/ + 3) 2 



[(2/ - l)P, +2 (o) - 2(2/ + l)P,(o) + (21 + 3)P_ 2 (o)] 



x [(2/ - l)Q l+2 (b) - 2(sl + 1)Q,(6) + (2/ + 3)Q ; _ 2 (6)] P,(z). (B24) 



56 



From the definitions (1B12j) and the identities (1B7j) . we obtain 

(l + 2)(l + 1)1(1- 1) 



(2Z-l)(2Z + l)(2Z + 3) 



(21 - l)P l+2 (a) 
-2(21 + l)P,(a) 



(2/ + 3)Pj_ 2 (a) , (B25) 



together with the corresponding equation for Qf(b). Using these relations, Eq. (1B24j) takes 
the form (lB17j) . as claimed. 



Finally, from Eqs. (IB16[) and flBlTj) . with a = — , b = — , x = cos# as before, we obtain 
the closed form expression f[9Tj) for ip^: 



tp4 



m AA 



1/2 

§ 2 



4M r 4 r 



M M 2 (t 2 + t 2 ) + t 2 t 2 — M 4 — 4tt M 2 cos + 2M 4 cos 2 9 1 



— m 



AA r 



3 A 3 /2 A 2 sin 2 g 

4 r r 4 i? 5 



(B26) 
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